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Abstract 

We systematically study the possible signals at LHC of new vector-like quarks 
mainly coupled to the third generation. We consider heavy quarks T, B, X, Y of 
charges 2/3, —1/3, 5/3 and —4/3, respectively, in SU(2)^ isosinglets Tl,r, -Bl,,r, 
or isodoublets (TE)^, (XT)l,r or {BY)l i r. Analyses based on a fast detector 
■ simulation are presented for twelve different final states containing one, two, three 

or four charged leptons in several invariant mass regions, also considering various 
b quark multiplicities. It is shown that with the combination of the different 
' channels the new quarks can be identified and their charged and neutral decays 

established. The comparison among final states also shows that the single lepton 
one offers the best discovery potential at LHC. For heavy quark masses of 500 
GeV, the 5a discovery luminosities range from 0.16 fb _1 for a (XT)l,r doublet 

dj \ to 1.9 fb _1 for a Bl % r singlet. 

! 

<N ■ 

1 Introduction 

*o : 
m . 

Hadron colliders offer an excellent place to look for new quarks, as the top quark 
discovery [1] and its recent observation in single production [2,3] at Tevatron evidence. 
In the near future, the operation of the Large Hadron Collider (LHC) will open a new, 
yet unexplored mass range for new physics searches, in particular for new quarks heavier 
than the top. Despite the demanding environment, new quark searches at LHC will be 
relatively clean because they can be produced in pairs through strong interactions with 
a large cross section and, being rather heavy, their signals can be easily distinguished 
from the large background from top pair production and W plus jets. 

Although the possibility of a fourth standard model (SM) sequential generation 
has not yet been excluded [4-6] and partial wave unitarity allows for fourth generation 
masses up to 1 TeV [7], new quarks heavier than the top are generally expected to be of 
vector-like nature if they exist. For example, extra-dimensional models with in the 
bulk [8-10] predict a tower of charge 2/3 isosinglets T^, of which the lightest one can 
be light and have sizeable mixing with the third generation [11]. More recently, (T B) l,r 
and (XT)l^r isodoublets of hypercharges 1/6, 7/6 coupling to the third generation 
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naturally emerge [12, 13] in warped models implementing a custodial symmetry to 
protect the Zbb coupling [14]. Charge —1/3 isosinglets B LR are predicted in grand 
unification theories based on E 6 , one of the most widely studied groups [15, 16], in 
which one such fermion per family appears in the 27 representation. Little Higgs 
models [17-19] also introduce a new Tl,r isosinglet partner of the top quark which 
ameliorates the quadratic divergences in the Higgs mass. In general, the new quarks 
predicted in these SM extensions are expected to couple mainly to the third generation. 
For generic Yukawa matrices and heavy quark mass terms, it has been shown [20] that 
the mixing of new vector-like quarks is of order m/M, where m,M are the masses of 
SM and new quarks, respectively. Then, unless specific symmetries are imposed on the 
mass matrices, the large mass hierarchy m t 3> m u ^, m& ^> jS favours mixing with 
the third generation. Additionally, constraints on top couplings are weaker than for 
the rest of quarks [21,22] so there is more room for mixing also from the experimental 
side. Note, however, that in some models it is possible to evade direct constraints with 
cancellations, and have large mixing with the first and second generations compatible 
with experimental data, see Ref. [23]. 

In case that any new physics is discovered at LHC, as it is hoped, it will be com- 
pulsory to determine its nature. For heavy quarks this means not only the observation 
of an event excess or even an invariant mass peak, but the determination of the quark 
charges and SU(2)^ isospin, the investigation of the decay channels and the measure- 
ment of their mixing with the SM quarks. In this paper we address some of these 
issues. We study the pair production of vector-like singlets T LjR , B LR of charges 2/3, 
—1/3 and doublets (TB)l,r, (XT)l,r, (BY)l^r of hypercharges 1/6, 7/6, —5/6, re- 
spectively, with quarks X, Y of charges 5/3, —4/3. (From now on we will drop the 
L, R subscripts.) We will assume that the new quarks mainly couple to the third 
generation. Previous literature has also investigated some of these signals in specific 
final states. For example, pair production of T singlets has been studied in the single 
lepton final state [24-26], as well as pair production of charge 5/3, —1/3 quarks in 



{XT), (T B) doublets producing like-sign dileptons [27] and one charged lepton [28] HI 
Here we will advance beyond previous work by analysing twelve multi-leptonic final 
states which give evidence for the several decay modes 



1 The discovery potential for D singlets coupling to u, d instead of the third generation has already 
been explored, for example in Refs. [29-31]. 




T W + b , T -> Zt, T -> Ht, 
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B -> Zb, B -> Hb, 

(1) 

with the aim of model discrimination. It is well known since some time (see for example 
Ref. [32]) that the presence or absence of specific decay modes can characterise the new 
quarks eventually observed. Here we demonstrate how this could be done in practice. 
For example, T quarks in a (X T) doublet have suppressed decay T — > W + b, so they 
are not seen in the W + bW~b final state as T singlets are. But they have enhanced 
T — > Ht decays, so if the Higgs boson is light (as preferred by electroweak precision 
data) they give a fairly large and clean TT — > Ht Hi — > HW + b HW~b signal with 
one charged lepton and six b quarks. On the other hand, YY — > W~bW + b cannot 
be distinguished from TT — > W + bW~b unless the b jet charge is measured, which is 
very difficult and requires large statistics. But, apart from different signal branching 
ratios, T quarks are cleanly identified by their characteristic T —> Zt decay, which can 
be observed in the trilepton final state. X and B quarks can both decay into four W 
final states, but in some models the latter also decays B — > Zb producing a sharp peak 
in a £ + £~b invariant mass distribution, which can be observed in dilepton and trilepton 
final states (here and in the following £ — e, /j,). In summary, here it will be shown that 
the simultaneous study and comparison of several multi-leptonic final states, with the 
observation of invariant mass peaks in most cases, can establish the identity of the new 
quarks, if they are observed at LHC. 

We remark that model discrimination is somewhat more demanding that evaluating 
the discovery potential of one's favourite model in some final state. From the technical 
point of view, it requires the complete signal generation with all decay channels. For 
TT and BB production there are in general nine decay modes according to Eqs. ([TJ). 
When the decay of the W and Z bosons (up to four, depending on the channel) are 
included, a plethora of possible final states appears involving multi-lepton signals. 
These contributions are all included in our simulations, which take into account the 
effects of radiation, pile-up and hadronisation, performed by a parton shower Monte 
Carlo, and use a fast detector simulation. SM backgrounds have also to be generated 
and simulated, including those with huge cross sections such as W and Z production 
plus jets, which are computationally demanding. 

Heavy quark pair production gives interesting signals in final states with one, two 
(like- and opposite-sign), three and four charged leptons. (Five and six lepton final 
states have too small branching ratios.) For model discrimination it is very convenient 
to classify signals not only by lepton multiplicity but by the number of Z boson "candi- 



B -> W~t, 
X -> W + t, 
Y -> W~b, 
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dates" present (same-flavour opposite-charge lepton pairs with an invariant mass consis- 
tent with M z ). For example, the trilepton final state is divided into a sample of events 
having a Z candidate (in which TT — > Zt W~b and other signals involving Z — > 
would be found) and events without Z candidates (to which XX — > W + t W~i, for in- 
stance, would contribute). In some cases the number of b jets present is also relevant. 
This gives a total of twelve interesting final states to be examined, and for which spe- 
cific analyses are presented in this paper. But, even after this final state organisation 
in terms of charged lepton multiplicity and number of Z candidates, there are final 
states where more than one type of quark give interesting signals. One of such cases 
is, for the trilepton final state with a Z candidate, 

TT -> Zt W~b -> Z W + b W~b Z -> tt , WW -> £vqq' , 

BB -> Zb W + t -> Zb W + W~b Z -> £+£', WW -»■ iuqq' (2) 

(the charge conjugate modes are also understood). In these likelihood classifi- 

cation is performed to separate and identify the TT and BB signals, and reconstruct 
them accordingly. This approach is unavoidable, since in some models like the [T B) 
doublet both signals can be present. 

Besides model discrimination, which is the main goal of this paper, the systematic 
study of all interesting final states offers several advantages. One of them is that the 
most sensitive ones can be identified. We find that the single lepton final state (with 
either two or four 6-tagged jets) offers the best discovery potential for all the models 
studied. For quark masses of 500 GeV, 5<r significance could be achieved for integrated 
luminosities ranging from 0.16 fb _1 for a (XT) doublet to 1.9 fb _1 for a B singlet. 
Our study also provides a guide of final state signatures to be searched in case that 
an event excess is identified in one of them. This complements previous work done 
for the characterisation and discrimination of seesaw models [33, 34] and new heavy 
leptons [35]. 

The structure of the rest of this paper is the following. In section [2] we introduce 
the models studied giving the relevant Lagrangian terms. In section [3] we discuss the 
general features of heavy quark pair production at LHC, and some details associated to 
the signal and background generation. In sections 0H8] the results for final states with 
four, three, two (like-sign and opposite-sign) and one charged lepton are presented, 
respectively. For the reader's convenience, the main results obtained are summarised 
at the end of each section, so that in a first reading the details can be omitted. Section [9] 
is a general summary where we address model discrimination by comparing signals in 
different final states. Our conclusions are drawn in section [TO] The Feynman rules 
used in our Monte Carlo programs are given in the Appendix. 
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2 Model overview 



In this section we briefly review the electroweak interactions of the new quarks, which 
determine their decay modes and single production. Additional details can be found 
in many early references, for example [20,32,36]. The interactions of (XT) and (BY) 
doublets are also given in Refs. [37,38]. 



2.1 T singlet 

We denote the SM weak eigenstates as q' u = (u' Li d' Li ) T ', u' Ri , d r Ri , where Latin indices 
i,j = 1,2,3 run over SM generations and Greek indices a, /3 = 1, ... ,4 over all quark 
fields. We use primes to distinguish them from mass eigenstates, where needed. The 
addition of a SU(2) L isosinglet u' Li , u' m to the SM quark content does not modify the 
SM charged and neutral current interactions in the weak eigenstate basis. The new u' m 
field has Yukawa couplings to the SM left-handed fields (the Yukawa coupling matrix 
Y must not be confused with a charge —4/3 quark Y), and a bare mass term can be 
written involving the new left-handed singlet u' L4 £ 

£z = [u' L iYu' Li - 2s 2 w J£ M ] , 

C Y = -Y$ + H.c. , 

Awe = -Mu' L4 u' m + H.C. (3) 

In this and the rest of models, the electromagnetic current J^ M has the same expression 
as in the SM but summing over all quark fields. The Higgs doublet is 

with v = 246 GeV and r the Pauli matrices. In the Lagrangians above we have omitted 
the terms in the down sector which are not affected by mixing. After the mass matrix 
diagonalisation the W, Z and H interactions read 

C w = -^u La 'fV a jdLjW+ + H.c, 

£z = -7T— [uLal^apULp ~ 2s w j£ M ] Z ^ , 

£v = ~WTT~ [uLaXapm u u R p + u Ra m u a X a pu Lp \ H, (5) 



2 In full generality, the right-handed fields u' Ra can be redefined so that the bare mass term only 
involves u' m . This change of basis also redefines the arbitrary matrix Y of Yukawa couplings. 



5 



where V a j is the 4x3 generalisation of the Cabibbo-Kobayashi-Maskawa [39,40] (CKM) 
matrix, X = VV^ a Hermitian 4x4 matrix (not to be confused with a charge 5/3 quark 
X) and the up-type quark masses. The electromagnetic current J^ M obviously re- 
mains diagonal. These equations, which result from a trivial change from weak to 
mass eigenstate basis, are exact and do not assume small mixing. Notice the appear- 
ance of left-handed flavour-changing neutral (FCN) couplings among up-type quarks, 
due to the mixing of left-handed weak eigenstates of different isospin, which breaks the 
Glashow-Iliopoulos-Maiani [41] mechanism. For a heavy quark T = 114 mixing with the 
top quark, and assuming small mixing, we have the approximate equality — Vm 
among neutral and charged current couplings, replacing generation indices by quark 
labels. This is a very well known result: in the T singlet model charged current mixing 
(WTb) automatically implies neutral current (ZTt) and scalar (HTt) interactions, all 
of the same strength up to multiplicative factors independent of mixing. The corre- 
sponding Feynman rules are given in the Appendix. These interactions determine the 
T quark decays, 

T^W + b, T^Zt, T —> Ht. (6) 
This new eigenstate has a mass tjit = M + 0(v 2 Y 2 /M 2 ). 



2.2 B singlet 

The Lagrangian for a B singlet is completely analogous to the one for a T singlet, with 
few replacements. The relevant interactions in the weak eigenstate basis read 

£z = -^[-d' Li rd' Li -2s 2 w J£ M ]Z lx , 
C Y = q' Li d' R p + H.c. 

Awe = -Md' LA d' m + H.C. (7) 

After mass matrix diagonalisation, we have 

Lw — 



■j= uu-fVipdLp W+ + H.c. , 



£>Z = —7T~ [—dLal^Xapdlfl — 1S 2 W J^ M ] Z ^ , 
ZCyy 

^ H = ~ oT7 [dLa^a/3 ™%d R p + d Ra m d a X a pd LI 3\ H . (8) 

ZMw 

The CKM matrix has dimension 3 x 4, X = V^V in this case and m d a are the down-type 
quark masses. For B mixing with the third generation we have X^ ~ VtB, so that the 
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new quark B has WtB, ZbB and HbB interactions governed by a single mixing factor 
VtB, m analogy with the T singlet model. The new quark B has a mass tub — M, and 
its decays are 

B -> W~t, B^Zb, B^Hb. (9) 



2.3 (TB) doublet 

With the addition of a vector-like doublet, the relevant Lagrangian in the weak inter- 
action basis is 

C w = -^[^7^L + ^47%4]< + H.c., 

£z = [u' La Yu' La + u' m Yu' m - d' La Yd' La - d' m Yd' m - 2s w J£ M } Z, , 

C Y = -Ylj q' La u' Rj <j> - Y d aj q' La d' Rj + H.c. , 
£ b are = -Mq LA q' m + H.c. , (10) 

with four SM-like left-handed doublets q' Li and one new right-handed doublet q' m = 
{u' m d' m ) T . The left-handed fields can be redefined so that the bare mass term only 
couples q' L4 . In the mass eigenstate basis it is more transparent to write the Lagrangians 
at first order in the (small) light-heavy mixing, 

C W = [uul^dLj + T L YB L + URa-fV^dRf,] W+ + H.c. , 

C Z = -■^^[uLa'fULa + URa'fX^URfi 

—dlal^dLa — dRal^ d a pdni3 — 2s w Jglvl] ' 

Ch = -^^[u L arnl{5 a p-X^)u Rp + u Ra {5 a p-X^)m u p u Ll3 

+d La m d a {5 aP - X d ap )d Rp + d Ra {5 aP - X d ^m d d Lp \ H , (11) 

so that it is apparent that the mixing of the heavy quarks T, B with SM quarks is 
only right-handed. The 4x4 matrix V R is not unitary, and also determines the FCN 
interactions, because X" = V^V^, X d = V RJf V R . Both X" and X d are Hermitian and 
non-diagonal, mediating FCN currents. Then, charged current interactions of the new 
states with SM quarks imply FCN ones, which result from the mixing of right-handed 
weak eigenstates with different isospin. At first order we have X tT ~ V^, X Bb ~ V R b . 
The new quarks are almost degenerate, with masses ttit = m R = M, up to terms 
of order v 2 Y 2 /M 2 . One can distinguish three scenarios for the heavy quark decays, 
depending on the relative sizes of the charged current mixing of the new quarks. For 
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V Tb ~ V t B the decay modes, assuming that they couple to the third generation, are the 
same as for singlets, 

T -> W + b , T —> Zt , T —> Ht, 

B^W't, B^Zb, B^Hb, (12) 

but with couplings of different chirality, which is reflected in some angular distributions. 
For Vrb "C VtB the top quark mixes with its partner much more than the bottom 
quark), the decays are 

T^Zt, T —> Ht, 

B -> W~t. (13) 

This scenario is the most natural one for generic Yukawa couplings due to the fact 
that the top quark is much heavier than the bottom quark, and is realised in some 
models [42]. Finally, a mixing Vrb VtB would give 

T -> W + b, 

B^Zb, B^Hb, (14) 

with signals similar to a hypercharge —5/6 doublet (BY) (see below). However, a 
mixing Vrb V t B is not natural in view of the mass hierarchy m t ^> m&, and is 
disfavoured by constraints on b quark mixing. 



2.4 (XT) doublet 

The interactions when a hypercharge 7/6 doublet is added have some similarities and 
differences with the previous case. In the weak eigenstate basis we have 
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JL 

£y = -Y« q' Li u' Rj cf> - Yl (X L u' L4 ) u' Rj <f> + H.c. , 

£bare = ~M (X L u' L4 ) I X * ) + H.C. , (15) 



£z = [u'txfu'u - u' L4 Yu' L4 - u' ml »u' m + X^X - 2s 2 w J£ M ] Z, , 



where for the charge 5/3 quark X the weak interaction and mass eigenstates coincide. 
We omit terms for the down sector which are unaffected by the presence of the new 
doublet. In the mass eigenstate basis, at first order in the light-heavy mixing the 
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Lagrangians read 

C w = [uuYVijdLj + X L YT L + XR-fVfpURp] W; + H.c. , 

C z = -7T- [uLiYu Li - f L -fT L - Una-fX^URf, + X^X - 2s 2 w J£ M ] , 

C h = -tJt~ [uL a m^(5 aP - X a/3 )u R/3 + u Ra (5 a p - X a/3 )mlu L p] H , (16) 

so that again the interactions of the new quarks X, T with the SM ones are right- 
handed. V L is the usual CKM matrix. The 1x4 matrix V R also determines the 
neutral mixing because X = V^V K . Notice an important difference with the T singlet 
and (T B) doublet: at first order the quark T does not have charged current couplings 
to SM quarks but has neutral ones ZTt, HTt. For a mixing with the third generation 
we have Xy t ~ V^ t . Obviously, the charge 5/3 quark X only has charged current 
interactions with SM charge 2/3 quarks. As in the (T B) doublet, the new mass 
eigenstates are almost degenerate, with masses m x — m T ~ M. Their allowed decays 
are 

X -> W + t, 

T^Zt, T —> Ht. (17) 



2.5 (BY) doublet 

Finally, the relevant Lagrangian for SM quarks plus a (B Y) doublet is 

C w = -^[u' L ^d' Ll + d' L ^Y L + d' m ^Y R ]w; + H.c. , 



£z = [-d' Li rd' Li + d' L4 Yd' L4 + d' m Yd' m - Y-fY - 2s w J£ M ] Z, , 

C Y = —Yfj q' Li d' Rj (f) — Yfj (d' L4 Y L ) d' Rj + H.c. , 

£ b are = ~M (d' L4 Y L ) ^ ^ + H.C. (18) 

At first order, the interactions in the mass eigenstate basis read 

C w = [u^YVfjdLj + B L YY L + dRa-fV^Yn] W+ + H.c. , 
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[-du-fdu + Bl^B l + d Ra ^X a/3 d R/3 - Y^Y - 2s 2 w J£ M ] Z^ , 



2cw 

C h = ~^h~ \dL a m d a (5 a p - X a/3 )d Rl3 + d Ra (5 a p - X a/3 )m d gd Ll3 ] H . (19) 

ZM w 

The matrix has dimension 4x1 and X = V R V R ^ . At first order the quark B does 
not have charged current couplings to SM quarks but has neutral ones. (The charge 
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—4/3 quark Y has only charged current interactions with down-type SM quarks.) For 
a mixing with the third generation we have X bB ~ Vf y . The new quarks have masses 
tub — my — M, and their allowed decays are 



B 
Y 



Zb, 
W~b. 



B -> Hb. 



(20) 



Notice that the Y — > W + b decay is like T — > W + b but with a b antiquark instead of a 
quark. These decays can be distinguished using angular distributions but, except for 
small kinematical differences, the signatures of a (B Y) doublet are similar to the ones 
of a (T B) doublet in which the B quark mixes much more than the T quark. 



3 Heavy quark production at LHC 



New heavy quarks can be produced in pairs via QCD interactions, 

gg,qq^QQ (Q = T,B,X,Y), 



(21) 



in the same way as the top quark. The cross section only depends on the quark mass, 
and is plotted in Fig. CD (left). For T quark singlets the partial decay widths are 



being 
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(23) 

a kinematical function. For a B singlet, the expressions for B — > VF~t, 5 — > Zb, B — > 
Hb can be obtained from Eqs. (l22l by replacing the mixings Vxb — > V t B, X Tt — > X Bb 
and the quark masses rriT — > ms, m t — > m b , m b — > m t . The branching ratios as a 
function of the heavy quark mass are presented in Fig. [T] (right), fixing Mh = 115 
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GeV. For (T B) doublets the analytical expressions of the widths are the same as for 
the singlets, although the relation beween the neutral and charged current mixings 
differs. For equal mixings Vrb — VtB the branching ratios are the same as for singlets, 
while for Vrb "C VtB the decays T — > W + b, B — > Zb, B — > Hb are absent, so that 
Br(T -> Zt) ~ Br(T -> if*) ~ 0.5, Bi(B -> Wi) = 1. For T, 5 quarks in (XT) and 
(£>F) doublets the charged decay modes are absent, and thus the partial widths for 
the other modes are roughly one half. For X — > W + t and Y — > W~b the widths are as 
for B — > W~t, T — ► replacing the mixings by and V^, respectively, as well 
as the quark masses. These are the only decay modes for X, Y quarks. 

Electroweak single heavy quark production is also possible at LHC, for example in 
the t-channel processes 

gq -> Tbq' , gq -> Tiq , 
gq -> fifog , -> Xtg' , 

gq^Ybq'. (24) 

For T6j, Xtj and F6j production (j = q,q' denotes a light jet) the processes involve 
a t-channel W boson, while Bbj and Tij production exchange a Z boson. This latter 
process has a much smaller cross section than Tbj but is the only possibility for the 
T quark in a (XT) doublet. The cross sections for the processes in Eqs. (f24l) are 
also plotted in Fig. [fl for reference mixings V,X = 0.1 with the third generation and 
including heavy quark and antiquark production. For the 2^2 processes bq — > Tj, 
bq — > Bj, etc. the cross sections are very close to the ones for their 2 — > 3 counterparts 
Tbj, Bbj, etc. Next-to-leading order corrections [44,45] are not included. In this work 
we do not consider single production as a means for model discrimination. Nevertheless, 
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depending on the heavy quark mass and mixing, some single production processes can 
be important, as it can be observed in Fig. Q3 In any case, it is important to remark 
that single production processes are crucial to measure the heavy quark mixing with 
SM quarks. 

The heavy quark signals studied in this paper have been calculated by implementing 
pair (Tf , BB, XX and YY) and single (Tj, Tbj, Ttj, Bj, Bbj, Xij, Yj, Ybj) 
production in the generator Protos [43], for the six models considered. All the decay 
channels in Eqs. (DQ) are included, with the subsequent W and Z boson decays in all 
channels. The Higgs boson decay, which does not carry any spin information, is left 
to the parton shower Monte Carlo. A complete signal evaluation is necessary for a 
study like the one presented here, which surveys final states from one to four leptons, 
and various b quark multiplicities in some cases. But also because sometimes charged 
leptons are missed by the detector, e.g. in Z — > resulting in contributions with 

fewer detected charged leptons than were generated at the partonic level. Matrix 
elements are calculated using HELAS [46], to take finite width and spin effects into 
account, and integration in phase space is done by Vegas [47]. The output is given in a 
suitable form to be interfaced to the parton shower Monte Carlo Pythia 6.4 [48] to add 
initial and final state radiation (ISR, FSR) and pile-up, and perform hadronisation. 

In this work we restrict our detailed simulations to heavy quark pair production, 
assuming heavy quark masses of 500 GeV and m t = 175 GeV, M H = 115 GeV. Cross 
sections and branching ratios are independent of the heavy-light mixing for T, B sin- 
glets and (XT), (BY) doublets, and a mixing V = 0.1 is assumed for definiteness. 
For the (T B) doublet we study two scenarios: (1) equal mixing Vrb = VtB = 0.1, in 
which the signals are quite similar to the ones of two T, B, singlets; (2) doublet mixing 
mainly with the top quark, Vrb = 0, V t B = 0.1. The signals for a doublet mixing 
mainly with the bottom are practically the same (except for the exchange of b quarks 
and antiquarks and small kinematical differences) as for the (BY) doublet, and are 
not presented for brevity. These six models are identified by the labels T s , B s , TB Al , 
TB d2 , XT A and BY A in tables and figures. Signals are generated with statistics of 300 
ftr 1 and rescaled to a reference luminosity of 30 fb -1 , in order to reduce statistical 
fluctuations. The factorisation and renormalisation scales used equal the heavy quark 
mass. We use the fast simulation AcerDET [49] which is a generic LHC detector simu- 
lation, neither of ATLAS nor of CMS, with standard settings. In particular, the lepton 
isolation criteria require a separation AR > 0.4 from other clusters and a maximum 
energy deposition HEt = 10 GeV in a cone of AR = 0.2 around the reconstructed 
electron or muon. Jets are reconstructed using a cone algorithm with AR = 0.4. In 
this analysis we only focus on central jets with pseudo-rapidity \r]\ < 2.5. Forward jets 
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with 2.5 < | T] | < 5 can also be present but are not considered for signal reconstruction 
nor for background rejection. For central jets, a simple b tagging is performed with 
probabilities of 60% for b jets, 10% for charm and 1% for light jets. We remark that 
the inclusion of radiation and hadronisation effects, as well as a detector simulation, is 
essential for our study. In an ideal situation in which the number of jets matches the 
number of partons in the hard process, the combinatorics to reconstruct the signals is 
relatively simple. In a real experiment, however, the presence of several more jets than 
were present at the partonic level, the radiation and the presence of mistags make it 
much more difficult to reconstruct and identify signals than it would be apparent with 
a toy parton-level simulation. An explicit example of these difficulties will be found in 
the single lepton channel in section [8], where we will show that TT and BB signals can 
sometimes be very alike, despite the very different decay chains involved. 

An adequate background calculation is another essential ingredient for our eval- 
uations. For multi-lepton signals, especially trileptons and like-sign dileptons, tinj 
(where nj stands for n additional jets at the partonic level) is one of the largest and 
most dangerous backgrounds, due to its large cross section and the fact that b quark 
decays sometimes produce isolated charged leptons. This background simply cannot 
be estimated with a parton-level calculation. Another important effect to be taken into 
account is the correct matching between the "soft" radiation generated by the parton 
shower Monte Carlo and the "hard" jets generated by the matrix element generator. In 
order to have predictions for SM backgrounds as accurate as possible we use Alpgen [50] 
to generate hard events which are interfaced to Pythia using the MLM prescription [51] 
to perform the matching avoiding double counting. The processes generated are col- 
lected in Table [Q where we also give the equivalent luminosity generated (30 fb _1 in 
most cases) and the number of events after matching. The additional SM processes 
bbnj and ccnj, which were previously shown to be negligible after selection cuts for 
multi-lepton states [33] are ignored in this work. (They are very likely to be negligible 
in the single lepton channel too, after the transverse energy and invariant mass cuts.) 

The procedure used for estimating the statistical significance of a signal is considered 
case by case. To claim discovery we require both (i) a statistical significance larger tan 
5a; (ii) at least 10 signal events. In the absence of any systematic uncertainty on 
the background, the statistical significance would be S Q = S/ y/B, where S and B 
are the number of signal and background events, or its analogous from the P-number 
for small backgrounds where Poisson statistics must be applied. Nevertheless, there 
are systematic uncertainties in the background evaluation from several sources: the 
theoretical calculation, parton distribution functions (PDFs), the collider luminosity, 
pile-up, ISR and FSR, etc. as well as some specific uncertainties related to the detector 
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Process 




Decay 


L 


Events 


tinj, n = 0, . . . 


,6 


semileptonic 


30 fb" 1 


6.1 M 


ttnj, n = 0, . . . 


,6 


dileptonic 


30 ft)" 1 


1.5 M 


tj 




->■ Zi/ 


30 ftj- 1 


0.9 M 


tb 




->■ Zi/ 


30 ft)" 1 


54 K 


tw 




all 


30 ft)" 1 


1.6 M 


titi 




all 


30 ft)" 1 


160 


tibb 




all 


30 ft" 1 


34 K 


Wnj, n = 0, 1, 


2 


Zi/ 


3 ft)- 1 


167 M 


Wnj, n = 3, . . 


.,6 


->■ Zi/ 


30 ft)- 1 


10 M 


Wbbnj, n = 0, 


...,4 


W 


30 ft" 1 


520 K 


Wccnj, n = 0, 


...,4 


W 


30 ft" 1 


550 K 


Wtinj, n = 0, . 


...,4 


W 


30 ft)- 1 


5.1 K 


n = 0, 1, 2, m« < 120 GeV 


z z+z- 


3 ft)- 1 


16.5 M 


Z/jnj, n = 3, 


. . . , 6, m H < 120 GeV 


z z+z- 


30 ftj- 1 


1.1 M 


Z/inj, n = 0, 


. . . , 6, m H > 120 GeV 


z -»• z+z- 


30 ft)- 1 


1.7 M 


Zbbnj, n = 0, . 


..,4 


z z+z- 


30 ft)- 1 


200 K 


Zccnj, n = 0, . 


..,4 


z -► z+z- 


30 ft" 1 


180 M 


Ztinj, n = 0, . 


..,4 


z i+i- 


30 fb- 1 


1.9 K 


WWnj, n = 0, 


...,3 


W 


30 ft)- 1 


290 K 


VFZnj, n = 0, 


...,3 


W ->lv,Z -> l + l" 


30 ft)- 1 


37.7 K 


ZZnj, n = 0, . 


..,3 


z z+z- 


30 ft)- 1 


3.7 K 


WWWnj, n = 


0,...,3 




30 ft)- 1 


1.5 K 


WWZnj, n = 


0,...,3 


all 


30 ft)- 1 


4.9 K 


WZZnj, n = 0,...,3 


all 


30 fb" 1 


1.5 K 



Table 1: Background processes considered in the simulations. The second column 
indicates the decay modes included (where / = e, /x, r), and the third column the 
luminosity equivalent generated. The last column corresponds to the number of events 
after matching, with K and M standing for 10 3 and 10 6 events, respectively. 

like the energy scale and b tagging efficiency. Such uncertainties have little relevance in 
the cleanest channels, where the discovery luminosity is controlled by the requirement 
of at least 10 signal events, being the significance far above 5a. For the channels 
in which the background normalisation can be important, we consider whether the 
signal manifests as a clear peak in a distribution. In such case it would be possible 
in principle to normalise the background directly from data, and extract the peak 
significance. Otherwise, we include a 20% background uncertainty in the significance 
summed in quadrature, using as estimator iS 2 o = S/ \JB + (0.2B) 2 . 
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4 Final state £+£+£-£- 



We begin our survey of the relevant final states with the one containing four leptons, 
which is the cleanest and less demanding one. The heavy quark signal reconstruction is 
not possible in most cases, but a simple event counting in several four-lepton subsamples 
already provides an useful test of the heavy quark signals. Having a small branching 
ratio in general, four leptons can be produced in several cascade decays of heavy quark 
pairs, for example 

Tf ^ ZtW'b ^ ZW + bW~b Z -> £ + £~,W -> Iv , 

Tf ^ZtVi^ ZW + b VW-b Z -> £ + £~, W -> tv, V -> qqjvv , 

BB^ZbZb Z^£ + £~, 

BB — > ZbW + i — > ZbW + W~b Z ^£ + £~,W ^£v, 

BB — ► W~tW + i — »■ W~W + bW + W~b W -> &/ , 

II -> W+tW-f -> W^->£z/, (25) 

with V = Z,H. The charge conjugate channels are implicitly included as well. The 
SM background is mainly constituted by ZZnj, tinj and Zttnj. The first one can be 
suppressed simply by requiring the presence of at least one 6-tagged jet, which hardly 
affects the signals which have two or more 6 quarks. Thus, for signal pre-selection 
we demand (i) four leptons summing a zero total charge, two of them with transverse 
momenta pt > 30 GeV and the other two with px > 10 GeV; (ii) at least one 6-tagged 
jet with pt > 20 GeV. We then develop three different analyses with disjoint event 
samples, aiming to separate the different signal sources of four leptons (ZZ, ZWW or 
WWWW leptonic decays). The criterion for the subdivision is the number of same- 
flavour, opposite-charge lepton pairs with an invariant mass consistent with Mz within 
some given interval, and the samples are labelled as 'ZZ\ 'Z J and 'no Z\ respectively. 
The invariant mass distribution of opposite-sign pairs can be studied by choosing pairs 
£^£^ as follows: 

1. If the charged leptons can be combined to form two Z candidates (there are two 
possibilities to construct two opposite-sign pairs), we label these pairs as t^l^ ■> 
£t^di ordered by transverse momentum. 

2. If not, we still look for a Z candidate combining opposite-sign pairs (there are four 
possible combinations). If found, we label this pair as £^£^ and the remaining 
leptons as £+, £j- 
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3. If no Z candidates can be found, we construct pairs £^£ b , £+£ d ordered by trans- 
verse momentum. 



The interval chosen to accept a Z boson candidate is M z ± 15 GeV, which provides a 
good balance between signal efficiency (for true Z boson decays) and rejection of non- 
resonant W + W~ decays giving opposite-charge leptons. The £^£7 and invariant 
mass distributions are presented in Fig. [2] for the six models, which are identified by 
the labels T s , B Sl TB dl , TB d2 (corresponding to the two mixing scenarios defined in 
section 12.31 ). XT d and BY d . These plots illustrate the relative size of the different 
signal contributions. Most signal events have at least one Z boson candidate: the 
WWWW decays correspond to £+£^ outside the Z peak (left plot). Events with two 
Z candidates are the ones with £^£ d at the Z peak (right plot). The distribution of 
signal and background events in the three samples at pre-selection is given in Table El 
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Figure 2: £+£ b , £~t£ d invariant mass distributions for the six models in the 
final state (see the text). The luminosity is 30 fb _1 . 





Total 


ZZ 


z 


no Z 




Total 


ZZ 


Z 


no Z 


Tf (T s ) 


50.0 


4.7 


33.3 


12.0 


BB (B s ) 


58.9 


12.3 


32.2 


14.4 


Tf (TB dl ) 


52.4 


3.9 


35.2 


13.3 


BB (TB dl ) 


54.3 


12.4 


28.3 


13.6 


Tf (TB A jXT d ) 


114.8 


12.1 


77.5 


25.2 


BB (TB d2 ) 


86.3 


1.2 


19.7 


65.4 


XX (XT d ) 


81.9 


1.0 


21.2 


59.7 


BB (BY d ) 


46.7 


29.7 


14.7 


2.3 












YY (BY d ) 


0.0 


0.0 


0.0 


0.0 


tinj 


7 





3 


4 


Ztinj 


15 





15 





Zbbnj 


1 





1 





ZZnj 


2 


2 









Table 2: Number of events in the £ + £ + £ £ final state for the signals and main back- 
grounds with a luminosity of 30 fb -1 , at pre-selection level. 
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4.1 Final state £+£+£-£- (ZZ) 

In this sample we do not impose any further requirement for event selection because 
the background is already tiny. The numbers of signal and background events can 
be read from Table El We observe that this final state is most useful for the model 
with a (BY) doublet where the decays B —> Zb are enhanced. The presence of the 
heavy quark B can be established by constructing a plot with the invariant mass of 
the 6-tagged jet and each of the two reconstructed Z bosons (two entries per event). 
This is shown in Fig. [3] for the six models considered, summing the contribution of 
the two quarks in the case of the doublets. The background, only two events, is not 
included. Notice that the bumps around 200 GeV in the TB d2 and XT d models cannot 
be mistaken by a charge —1/3 quark even with low statistics: for such a mass the heavy 
quark production cross section would be more than 100 times larger. 




Figure 3: £^£ b b, £+£ d b invariant mass distribution for the six models in the £ + £ + £ £ 
(ZZ) final state, with two entries per event. The luminosity is 30 fb _1 . 

We give in Table[3]the luminosity required to have a ha discovery, including all signal 
contributions within a specific model. The background normalisation uncertainty has 
little relevance in these cases, because the background itself is very small and the 
discovery luminosity is mainly determined by the minimum of 10 signal events. We 
also include in this table whether a mass peak can be reconstructed, although in this 
case the peak observation and reconstruction clearly requires more luminosity than ha 
discovery, due to the small statistics. We point out that, since the heavy resonance 
is observed to decay into a Z boson and a b quark, it can be identified as a heavy B 
quark. This, however, can also be done in the opposite-sign dilepton final state with 
six times better statistics. For the models with T quarks, mass peaks could in principle 
be reconstructed with a high integrated luminosity, but this is far more interesting to 
do in the trilepton channel where statistics are larger. 
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L Rec. 

T s no 

B s 24 ft^ 1 m B 

TB dl 18 fb" 1 m B 



L Rec. 

TB d2 23 fb" 1 no 

XT d 23 fb" 1 no 

BY A 10 fb" 1 m B 



Table 3: Luminosity L required to have a 5a discovery in the £ + £ + £~£~ (ZZ) final 
state. A dash indicates no signal or a luminosity larger than 100 fb -1 . We also indicate 
whether a mass peak can be reconstructed in this final state. 

4.2 Final state £+£+£-£- (Z) 

Events with only one Z boson candidate are selected in this sample. Additionally, 
the presence of two (6-tagged or not) extra jets with pt > 20 GeV is required to 
reduce the Ztinj background, hardly affecting the signals. The number of signal and 
background events at pre-selection and selection is collected in Table BJ Notice that 
for XX production, where Z bosons are not produced in the decay, in some cases a 
pair of charged leptons from W + W~ decays accidentally have an invariant mass in 
the interval selected. Nevertheless, this non-resonant contribution is 5 times smaller 
than the one from pair production of its T partner. The same comment applies to BB 
production in the TB d2 model. 





Pre. 


Sel. 




Pre. 


Sel. 


Tf (T s ) 


33.3 


29.5 


BB (B s ) 


32.2 


25.1 


Tf (TB dl ) 


35.2 


31.5 


BB (TB dl ) 


28.3 


21.4 


Tf (TB d2 /XT d ) 


77.5 


74.9 


BB (TB d2 ) 


19.7 


14.0 


XX (XT d ) 


21.2 


15.5 


BB (BY d ) 


14.7 


12.5 








YY (BY d ) 


0.0 


0.0 


tinj 


3 





Ztinj 


15 


8 


Zbbnj 


1 





ZZnj 









Table 4: Number of events in the £ + £ + £ £ (Z) sample for the signals and main 
backgrounds with a luminosity of 30 fb -1 . 

The reconstruction in this final state is very difficult due to the presence of two 
final state neutrinos, each resulting from the decay of one heavy quark. Thus, we 
restrict our analysis of this sample to the Z boson identification and a simple counting 
of events, which can already be an useful test of the different models. In Table [5] we 
collect the discovery luminosities for the six models studied. We observe that those 
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L 


Rec. 




L 


Rec 


Z 


11 fb" 1 


no 


TB d2 


3.4 fb" 1 


no 


B s 


14 fb" 1 


no 


XT A 


3.3 fb- 1 


no 


TB dl 


5.7 fb" 1 


no 


BY d 


50 fb" 1 


no 



Table 5: Luminosity L required to have a 5a discovery in the £ + £ £ + £ (Z) final state. 
We also indicate whether a mass peak can be reconstructed in this final state. 

with T quarks give important signals, especially the ones with enhanced branching ratio 
for T — > Zt, and the discovery luminosities are relatively small. In these interesting 
cases the background normalisation uncertainty is not important because the signals 
are much larger. 

4.3 Final state £+£+£-£- (no Z) 

This sample contains the signal and background events for which all opposite-sign pairs 

have invariant masses \m,+,- — M z \ > 15 GeV. We do not apply any further event 

* i 

selection criteria since the background at pre-selection is already rather small. The 
number of signal and background events can be read in Table El The most important 
signals are from XX production, for which the decay X — > W + t — > W + W + b has 
branching ratio unity, and BB production in the TB d2 model, with unit branching 
ratio for B — > W~t — > W~W + b. The latter decay approximately has a branching 
ratio of 0.25 for the B singlet and TB dl doublet models, and is absent for the (BY) 
doublet. TT production, which in the four-lepton final state at least involves one Z 
leptonic decay, gives a small contribution which is only due to the finite Z width and 
energy resolution of the detector. 





L 


Rec. 




L 


Rec 


T s 


35 fb" 1 


no 


TB d2 


3.3 fb" 1 


no 


B s 


25 fb" 1 


no 


XT d 


3.5 fb" 1 


no 


TB dl 


11 fb" 1 


no 


BY d 




no 



Table 6: Luminosity L required to have a ha discovery in the £ + £ + £~£~ (no Z) final 
state. A dash indicates no signal or a luminosity larger than 100 fb -1 . We also indicate 
whether a mass peak can be reconstructed in this final state. 

We collect in Table [6] the luminosity required for 5a discovery of the six models 
considered in this work. The reconstruction in this final state is virtually impossible 
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because four neutrinos are present in the final state and, in fact, all like-sign and 
opposite-sign dilepton distributions seem very similar. Nevertheless, as in the previous 
sample, the number of events itself is a very good check of the different models. For the 
most interesting signals (with (T B) and (X T) doublets) the background normalisation 
is not important, while for the other cases the luminosities given are a little optimistic. 

4.4 Summary 

Four lepton final states have seldom been considered in the context of heavy quark 
searches, perhaps because they are less relevant for the traditionally most popular 
models with T or B singlets. Nevertheless, for the (X T) and (B Y) doublets and 
the TB& 2 model the multi-lepton signals are larger in general: either for the decays 
T — > Zt, B — > Zb (which have branching ratios two times larger than in the singlet 
case), or from the decays X — > W + t, B — > W~t (with unit branching ratio). Thus, 
the four-lepton final state can be interesting for this class of models. One has to note 
here that the sensitivity to heavy quark signals in other final states is much better, and 
discovery luminosities one order of magnitude smaller. Still, four lepton signals would 
be visible with a moderate luminosity and should be explored to test the models. 

It is very convenient to divide the four lepton final state in three different subsets 
( L ZZ\ 'Z' and 'no Z y ) depending on the number of Z boson candidates (2, 1 and 
0, respectively) present. This subdivision allows for some model discrimination from 
event counting in this final state alone, for example: 

• If a signal is simultaneously observed in the 'Z' and 'no Z J samples with a similar 
luminosity, but not in the 'ZZ J one, it points towards a (XT) doublet or a (T B) 
doublet predominantly mixing with the top quark (TB d2 model). 

• If, conversely, a signal is observed exclusively in the 'ZZ J sample, it corresponds 
to a (BY) doublet. The presence of the heavy B quark can also be established 
by the observation of a peak in the Zb invariant mass distribution. However, this 
can also be done in the opposite-sign dilepton final state with six times better 
statistics. 

Finally, it is worth mentioning that the four lepton final state is also a possible signal 
of heavy charged lepton in several models [35], but in that case the invariant mass of 
three charged leptons displays a very clear and sharp peak at the heavy charged lepton 
mass rriE, and b quarks are not produced. Four leptons are also produced in the decay 
of doubly charged scalars produced in pairs (for a detailed analysis see Ref. [33]) but 
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for the scalar triplet signals are clearly distinguishable by the presence of narrow peaks 
in the like-sign dilepton invariant mass distributions. 

5 Final state PlH* 

The trilepton final state offers a good balance between signal branching ratio in TT, 
BB and XX production, and SM background. Three leptons can result from several 
heavy quark pair cascade decays, either involving the leptonic decay of a Z and a W 
boson, as for example in 

TT -> Zt W~b -> ZW+bW'b 
Tf -> ZtVi^y ZW + b VW'b 
BB -> ZbW + t -> ZbW + W'b 

with V — Z,H, or of three W bosons, 

BB -> W"t W + t -> "6 3iy -»■ 1W -> gg , 

XX -> lf + ir< -> W + W + 6W"W-6 — > £z/, 1W — > gg' . (27) 

The charge conjugate channels are implicitly included in all cases. All these production 
and decay channels are interesting and a first signal discrimination can be made, as in 
the previous section, by the presence or not of Z boson candidates in the final state. 
In the sample with Z candidates it is necessary to go further and try to separate the 
three channels in Eqs. (l26l) . An obvious reason motivating this separation is that for 
(T B) doublets both TT and BB pairs can be produced and the three processes in 
Eqs. (l26l) are present in general. Then, it is quite desirable to separate the signals of T 
and B quarks, identifying their production and decay channels. The discrimination is 
possible with a probabilistic analysis which classifies the events into the three processes 
in Eqs. (|26l) with a good efficiency. 

The main SM backgrounds to trilepton signals are from WZnj and itnj production, 
both roughly of the same size. The latter is originated when the two W bosons decay 
leptonically and one b quark gives a third isolated lepton but, as in the like-sign dilepton 
final state examined in the next section, it can be significantly reduced by asking that 
the two like-sign leptons have high transverse momenta. Thus, for event pre-selection 
we require the presence of three charged leptons (summing a total charge ±1), the like- 
sign pair having px > 30 GeV and the third lepton with pt > 10 GeV. As mentioned 
above, we divide the trilepton sample into two disjoint ones. The first one contains 
events where a Z boson candidate can be identified, that is, when two same-flavour 



z - 


-> £ + £~ 


,ww 


z - 


*£+£- 


,ww 


z - 


-> £ + r 


,ww 



(26) 
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opposite-charge leptons have an invariant mass consistent with Mz- The other sample 
contains events without Z candidates. The interval in which a lepton pair is accepted 
as a Z candidate is chosen to be of 15 GeV around Mz- We can compare the signal 
contributions to the two samples by plotting the invariant mass of two opposite-charge 
leptons £+, £^, chosen in the following way: 

1. If there is a Z candidate, we label the corresponding leptons as £ b . In case 
that there are two Z candidates, which can accidentally happen, the leptons with 
largest transverse momenta are chosen. 

2. If there are no Z candidates, we choose £^ with the largest transverse mo- 
menta. 




Figure 4: £+£ b invariant mass distributions for the six models in the tr^l^l^ final state 
(see the text for the definition of £+ and £ b ~). The luminosity is 30 fb _1 . 



The resulting distribution is shown in Fig. 01 We observe that there is a large off- 
peak signal from BB and XX decays in the TB d2 and XT d models, respectively. The 
number of events at pre-selection in each sample is given in Table [7J A sizeable fraction 
of events from TT decays, which in the trilepton channel always involve a Z boson, are 
classified in the 'no Z 7 set, while around 10% of the BB and XX events in which Z 
bosons are not present are accepted in the 'Z' sample. The rate of wrong assignments 
can be reduced at the cost of losing signal efficiency, by strenghtening the classification 
criteria. For example, Z candidates could be accepted only in the interval M^±10 GeV 
and events in the 'no Z' subsample could be rejected if opposite-charge pairs have an 
invariant mass in the range Mz ± 20 GeV. This fine tuning of the analysis makes more 
sense with a full detector simulation, and is not necessary for model discrimination, 
anyway. 
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Total 


Z 


no Z 




Total 


Z 


no Z 


Tf (T s ) 


320.7 


212.4 


108.3 


BB (B s ) 


421.9 


227.9 


194.0 


Tf (TB dl ) 


349.0 


229.9 


119.1 


BB (TB dl ) 


484.5 


237.0 


247.5 


Tf (TB d jXT d ) 


654.6 


435.8 


218.8 


BB (TB d2 ) 


1174.4 


144.0 


1030.4 


XX (XT d ) 


1181.8 


143.9 


1037.9 


BB (BY d ) 


106.3 


88.3 


18.0 










YY (BY d ) 


0.5 


0.1 


0.4 


tinj 


464 


114 


350 


WZnj 


4258 


4196 


62 


Wtinj 


78 


11 


67 


ZZnj 


424 


417 


7 


Ztinj 


189 


169 


20 











Table 7: Number of events in the t^t^t^ final state for the signals and main back- 
grounds with a luminosity of 30 fb -1 , at pre-selection level. 

5.1 Final state t^i^l^ (Z) 

This final state receives important contributions from Tf and BB production in the 
channels of Eqs. ([261 . We will first perform an analysis with fewer selection criteria 
to suppress the background and obtain the heavy quark discovery potential for this 
final state. Then, we will address the identification of a heavy quark signal eventually 
observed, strengtheining our requirements on signal and background events and using 
a likelihood function which assigns them to each of the decay channels in Eqs. fl26l) . 
After that, events will be reconstructed accordingly to their classification and to the 
kinematics assumed in each case. 

5.1.1 Discovery potential 

For events with one Z candidate we ask (i) at least two light jets with pt > 20 GeV; (ii) 
one 6-tagged jet also with pt > 20 GeV; (iii) transverse momentum pr > 50 GeV for 
the leading charged lepton l\\ (iv) transverse energy H T > 500 GeV. The kinematical 
distributions of these variables at pre-selection are presented in Fig. [5] for the relevant 
signals and the SM background. In particular, requiring a 6-tagged jet hardly affects 
the signals but practically eliminates the WZnj background which does not have b 
quarks. The cuts on transverse energy and leading charged lepton momentum are 
quite general to look for new heavy quarks and are not optimised for the input masses 
used in our calculation. Notice also that the Ht distribution for the signals clearly 
indicates that one or more heavy particles with masses summing around 1 TeV are 
produced. This data will be crucial later when we address the disentanglement and 
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reconstruction of different signal channels. 

The number of signal and background events at the selection level is given in Ta- 
ble El also including the values at pre-selection for better comparison. As it might be 
expected, the most important background after cuts is Ztinj, which has a Z boson, 
two b quarks, two W bosons and large transverse energy, and is then quite similar to 
the signals. More aggresive cuts will of course reduce this and the other backgrounds 
but we refrain ourselves from performing such optimisations. The discovery luminosi- 
ties are given in Table [9J summing all signal contributions within a given model. We 
observe that this clean channel offers a good potential to discover T and B quarks 
in singlet or doublet representations. For the (BY) doublet the discovery luminosity 
may be optimistic because it does not take into account the background normalisation 
uncertainty, which may be important in this case where the signal is small. 

5.1.2 Heavy quark reconstruction 

The broad sensitivity to TT and BB signals of this final state implies that if a positive 
excess is observed, identifying its nature will require a more elaborate analysis. Indeed, 
the decay modes in Eqs. ( 1261 ) give signals only differing by the number of jets and the 
location of the resonant peaks. The identification can be done efficiently, however, by 
using a likelihood method which gives the probability that a given event corresponds to 
each of the decay modes. We build probability distribution functions (p.d.f.) for three 
signal classes: (a) TT — > ZtWb; (b) TT ZtVt, with V = H, Z decaying hadronically 
or invisibly; (c) BB — > ZbWt. We generate high-statistics samples different from the 
ones used for the final analysis. We do not include a separate class for the background, 
because the likelihood function is only used to identify signals and not to reject the 
background which is rather small. Nevertheless, the discriminant analysis also affects 
the background and, for instance, if we restrict ourselves to events classified as result- 
ing from BB production, a sizeable part of the background is classified as TT-like and 
thus rejected. Note also that an essential parameter for building the kinematical dis- 
tributions for the signals is the heavy quark mass. If a heavy quark signal is observed 
at LHC, the approximate value of the heavy quark mass can be estimated from the 
transverse energy distribution for the signaljj and then a probabilistic analysis can be 
performed to separate signal contributions and reconstruct the decay chain event by 
event. 

For the signal discrimination and reconstruction we demand, in addition to the 

3 In the case of the doublets the new states are expected to be nearly degenerate, simplifying their 
approximate determination from this distribution. 
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Figure 5: Kinematical distributions of variables used in selection and recontruction 
criteria for the i^i^i^ (Z) final state: light jet multiplicity Nj, b jet multiplicity N b , 
transverse momentum of the leading lepton and total transverse energy. The luminosity 
is 30 ftr 1 . 
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Table 8: Number of events at the pre-selection, selection and reconstruction levels in 
the £ ± £ ± £ =F (Z) sample for the signals and main backgrounds with a luminosity of 30 
ft)" 1 . 
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Table 9: Luminosity L required to have a 5a discovery in the £ ± £ ± £ T (Z) final state. 
We also indicate whether a mass peak can be reconstructed in this final state. 

selection criteria already specified, the presence of at least two 6-tagged jets. The 
number of signal and background events with this last requirement is given in Table El 
We begin by finding two W bosons decaying hadronically and leptonically. The former 
is approximately reconstructed at this stage by selecting among the light jets with 
largest pt (up to a maximum of four) the two ones which give an invariant mass closest 
to Mw- The latter is approximately reconstructed from the charged lepton not included 
in the Z candidate and the missing energy, with the procedure explained below, and 
selecting the solution giving the smallest neutrino energy. The variables used in the 
likelihood function are: 

• The light jet and b jet multiplicities. 

• The invariant mass of the reconstructed W boson (decaying hadronically or lep- 
tonically) with larger transverse momentum, labelled as W\, plus the b quark 
with largest transverse momentum, b\. For T — > Wb decays the W boson as well 
as the b quark are expected to have larger pt, and we observe in Fig. [6] that this 
is often the case. 

• The invariant mass of the reconstructed Z boson and the b quark with highest 
Pt, which for the BB signal is most times the one resulting from B — > Zb, as we 
observe in the distribution of Fig. [6l 

• The invariant mass of the reconstructed W with smaller pr (W 2 ), the Z boson 
and the b quark with smaller transverse momentum (62)- 

• The invariant mass of the two W bosons and the b quark with smallest px, which 
for the BB signal are the ones from B — > Wt — > WWb in most cases. 

The likelihood function evaluated on the three class samples gives the probability dis- 
tributions in Fig. [71 where P a , P&, P c are the probabilities that events correspond to 
each of the three likelihood classes in Eq. (1261) . Events are assigned to the class (a, b 
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Figure 6: Kinematical variables used to classify the three heavy quark signals in the 
(Z) final state. 



or c) which has the highest probability P a , Pb or P c , respectively. Table fTOl shows the 
performance of the likelihood function on the reference samples. Events in a class x are 
correctly classified if P x > P y ,P z , where y, z are the other classes. The probabilities 
for correct assignments are in the range 0.61 — 0.69, which suffice to achieve a good 
reconstruction of the heavy resonances. We now describe the procedure followed in 
each case. 

Class (a): TT — > ZtWb —> ZWbWb. Events which are identified as resulting from 
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Figure 7: Probability distribution functions for events in the reference samples. 



Class P a >P b ,P c P b >P a ,P c P c >P a: P b 



(a) 


0.61 


0.24 


0.15 


(b) 


0.19 


0.69 


0.12 


(c) 


0.15 


0.20 


0.65 



Table 10: Performance of the likelihood function on the £ ± £ ± £ =F event reference samples: 
fractions of events in each sample and their classification. Events in a class x are 
correctly classified if P x > P y , P z , where y, z are the other classes. 

this decay chain are reconstructed using this procedure: 

1. The Z boson momentum is obtained from the opposite-sign lepton pair Z candi- 
date. 

2. Two light jets are selected to form the hadronic W, labelled as Wh- If there are 
only two light jets these are automatically chosen; if there are more than two, 
only up to four (ordered by decreasing p?) are considered. 

3. The leptonic W (labelled as Wl) is obtained from the charged lepton t not 
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included in the Z candidate and the missing energy, identifying (p u )r = v/r , 
requiring (p£ + p u ) 2 = and solving for the longitudinal component of the 
neutrino momentum. If no real solution exists, the neutrino transverse momen- 
tum is decreased in steps of 1% and the procedure is repeated. If no solution is 
still found after 100 iterations, the discriminant of the quadratic equation is set 
to zero. Both solutions for the neutrino momentum are kept, and the one giving 
best reconstructed masses is selected. 

4. Two b jets are selected among the ones present, to be paired with Wh and Wl, 
respectively. 

5. The top quark is reconstructed from one of the Wb pairs, and its parent heavy 
quark 7\ from the top quark and the Z boson. 

6. The other heavy quark T 2 is reconstructed from the remaining Wb pair. 

7. Among all choices for b and light jets and all possible pairings, the combination 
minimising the quantity 

(mjfc - M W f (mj% - M W ) 2 ( mf c _ m) 2 (m ™ _ m rec )2 

V 2 + V 2 + V 2 + V 2 — { ' 

is selected, with a w = 10 GeV, a t = 14 GeV [52], a T = 20 GeV. Notice that 
we include the leptonic W boson reconstructed mass in the minimisation. Since 
the quadratic equation is forced to have a solution in all cases, sometimes the 
reconstructed mass is not the W mass. 

Class (b): TT — > ZtVi — > ZWbVWb. For events identified as resulting from this 
decay chain which have at least six jets (otherwise they they are dropped) we proceed 
through the same steps 1 — 4 as in class (a), and subsequently: 

5. The hadronic and leptonic tops tn, ti are obtained from the two Wb pairs. 

6. One heavy quark Ti is reconstructed from one top and the Z boson. The other 
heavy quark is obtained from the other top and two jets chosen among the ones 
present (6-tagged or not). 

7. The combination minimising 

(m™ H - Mw ) 2 | (m™ - M W ) 2 | (m^ - m t ) 2 | {m\f - m t ) 2 | (mgf - m^) 2 

is finally selected. 
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Class (c): BB — > ZbWt — > ZbWWb. The reconstruction of this channel proceeds 
through the same steps 1 — 3 as in the previous two channels, and then: 

4. Two b jets are selected among the ones present, and one of them is paired with 
the Z boson to reconstruct a heavy quark B\. 

5. The second b jet is associated with one of the W bosons to reconstruct a top 
quark, and then with the other W boson to reconstruct the second heavy quark 
B 2 . 

6. The combination minimising 

(m^ H -M w f (m^ L -M w f | [m^ - m t f (m% - mg) 2 

„2 + „2 + „2 + „2 \ 6y> ) 

a W W t B 

is finally selected, with <tb = 20 GeV. 
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82.9 
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65.3 


15.3 
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41.7 


11.4 
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41.0 
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24.2 
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34.2 


12.3 


11.8 


3.2 


XX (XT d ) 


35.7 


14.2 


12.6 
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21.2 
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2.7 


11.3 












YY (BY d ) 


0.0 


0.0 


0.0 


0.0 



Zttnj 32 5 12 5 



Table 11: Number of signal and background events in the £ ± £ ± £ T (Z) final state at the 
reconstruction level assigned to each event class. The luminosity is 30 fb -1 . 

We present our results in Fig. El including all signal contributions in a given model, 
as well as the SM background, and discuss them in turn. We do not include W and 
top reconstructed masses, which show good peaks at the true masses with the opti- 
mised method used. The separate contributions of each process are given in Table [EH 
skipping several backgrounds which are practically removed at the last stage of event 
selection. (The total number of events includes in each case those in class (b) which 
are later rejected by the reconstruction algorithm because they do not have at least 6 
jets.) The first plot (up, left) shows the reconstructed T\ mass for events assigned to 
classes (a, b). This heavy quark is the one decaying T — » Zt, with t decaying either 
hadronically or semileptonically. The reconstruction of this peak in the Zt invariant 
mass distribution implies that T has charge 2/3, and also shows the vector-like nature 
of T. The counterpart for B quarks is shown in the second plot (up, right), with the 
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Figure 8: Reconstructed heavy quark masses in the i^i^i^ (Z) final state. 

reconstructed mass of B 1 , which is the quark decaying B — > Zb. The reconstruction 
of a peak in the Zb invariant mass distribution shows that B has charge —1/3 and is 
vector-like. 

The other resonant peaks also give information regarding heavy quark decays. We 
show in the third plot (middle, left) the reconstructed T 2 mass for events in class 
(a), which corresponds to the decay T — > Wb, with W decaying either leptonically or 
hadronically. For the T singlet and (T B) doublet in scenario 1, where this decay takes 
place, the peaks are sharp, and they might be observed with sufficient luminosity. We 
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point out that the presence of events with W decaying leptonically (about one half 
of the total) clearly indicate the T — > Wb decay, but this can also be established in 
the single lepton final state with much larger statistics. The other models with T 
quarks in which the decay T — > Wb does not happen still have a fraction of events 
incorrectly assigned to this class. In these models the Wb invariant mass distribution, 
which should peak at m tl is broader and shifted towards larger values because the 
reconstruction procedure enforces equal masses for both heavy quarks. The fourth plot 
(middle, right) represents the B 2 invariant mass distribution for events in class (c), 
from the decay B — > Wt — > WWb. This plot also shows the presence of a resonance 
decaying into a top quark and a W boson, the latter reconstructed either from two 
jets or from a charged lepton plus missing energy. This peak establishes the decay 
B — ► Wt, which is absent in the (B Y) doublet. Finally, the fifth plot (down) shows 
the reconstructed T 2 mass in class (b), corresponding to the decay T — > Vt, with V 
decaying into two jets. This distribution shows the presence of a resonance but does 
not help establish its nature, because the identity of V is not determined. 

A few remarks are in order. It is clear that detecting the presence of a resonant 
peak and drawing conclusions about the nature of the heavy quark requires a significant 
amount of statistics, and a compromise should be taken between having good recon- 
structed peaks (imposing quality cuts on class identification as well as on reconstructed 
W boson and top quark masses, for example) and having a sufficient number of events. 
Here we have made no quality cuts in order to keep the signals as large as possible. But 
even with this conservative approach the contributions of the three cascade decays in 
Eqs. ( l26l) can be disentangled, and invariant mass peaks can be reconstructed so that, 
if sufficient luminosity is collected, the decays T — > Zt, T — > Wb, B — ► Zb, B — > Wt 
can be established. 

Finally, we address the discrimination of T singlets and T quarks of a (T B) doublet 
in scenario 1, using angular distributions. (In this scenario the T decay branching 
ratios are the same as for T singlets.) In T — ► Zt decays, as well as in B —> W~t, 
the top quarks are produced with a high polarisation P = ±0.91 in the helicity axis 
(negative for the singlets and positive for the doublets), and the opposite polarisation 
for antiquarks. This allows to determine the chirality of the WTb and WtB couplings 
by looking at the charged lepton distribution in the top quark rest frame for the event 
subset in which the top decays leptonically. We show in Fig. [9] (left) the theoretical 
distributions as computed with the Monte Carlo generator, which are the same for T 
and B quarks, since the decays T — > Zt and B — ► W~t involve a coupling with the same 
chirality, left-handed for singlets and right-handed for doublets. On the right panel we 
show the reconstructed distribution for the T singlet and TB Al model, including in 
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Figure 9: Left: Charged lepton distribution in the top quark rest frame for T — > Zt 
and B — > W~t decays. Right: distribution for the T singlet and (T B) doublet after 
simulation. 

the latter case the B contribution which is flat and slightly smooths the slope of the 
distribution. It is clear that large statistics are required to discriminate both cases, 
but the differences are visible already without the need of unfolding detector effects. 
The forward-backward asymmetries computed from the reconstructed distributions are 
A F b = —0.19 for T s and A FB = 0.24 for TS dl , so that with 30 fb -1 (corresponding to 
27.7 and 38.6 events in each case) the statistical difference would amount to 2.4a. A 
complete analysis unfolding the detector effects and with an appropriate calculation of 
systematic uncertainties (see for example Ref. [53] for a similar analysis) is beyond the 
scope of this work. For B — > W~t decays the results are completely analogous but with 
smaller statistics. We also note that for these large heavy quark masses the Z bosons 
produced in T — > Zt, B — > Zb are mostly longitudinal, and the angular distribution of 
the pair from Z decay is almost indistinguishable for T, B singlets and doublets 
already at the generator level. 

5.2 Final state (no Z) 

In the sample without Z candidates we ask for event selection that (i) the leading 
like-sign lepton l\ has transverse momentum > 50 GeV; (ii) the total transverse 
energy H T is larger than 500 GeV. Notice again that these cuts are not optimised 
to reduce the background but are quite general to search for new heavy states. The 
kinematical distributions of the two variables at pre-selection are shown in Fig. [TO] for 
the SM background and all models except the (B Y) doublet model which has a very 
small signal. As in other final states, the H T distribution clearly indicates in all cases 
that one or more heavy particles, summing a mass around 1 TeV, are produced. The 
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Figure 10: Left: transverse momentum distribution of the leading like-sign lepton. 
Right: total transverse energy. The luminosity is 30 fb _1 . 
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Table 12: Number of events at the pre-selection, selection and reconstruction levels in 
the i^l^l^ (no Z) final state for the signals and main backgrounds with a luminosity 
of 30 fb" 1 . 

number of events after this selection is given in Table CE21 including also the numbers at 
pre-selection for comparison. We do not require 6-tagged jets at this stage because it 
does not improve the background rejection, since most of the backgrounds have two top 
quarks. The 5cr discovery luminosities in this final state for the six models, including 
the contribution of both members in the case of the doublets, are given in Table [131 
The luminosity indicated for the T singlet may be optimistic because the background 
normalisation uncertainty (unimportant for the other models which have much larger 
signals) may be relevant. 

In this final state the heavy quark masses cannot be directly reconstructed because 
each heavy quark has among its decay products an invisible neutrino, and there are 
three neutrinos in total. In XX production we can still have information about the 
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Table 13: Luminosity L required to have a 5a discovery in the £ ± £ ± £ =F (no Z) final 
state. A dash indicates no signal or a luminosity larger than 100 fb -1 . We also indicate 
whether a mass peak can be reconstructed in this final state. 

mass of one of the heavy quarks, the one decaying X — > Wt — > WWb, with WW — > 
iuqq', by reconstructing its decay products except the missing neutrinoo For the 
reconstruction we demand the presence of at least one 6-tagged jet and two light jets 
with px > 20 GeV, and select the hadronic X decay products as follows: 

1. We select a b jet among the ones present. 

2. We select a pair of light jets ji, j 2 among the three ones with highest px (in case 
there are only two, we select these. 

3. We choose the combination minimising the quantity 

(m jlj2 - M w ) 2 (m jlj2b - m t f 

2 + 2 (61) 

°w °t 

with a w = 10 GeV, a t = 14 GeV. 

The "visible" component of the heavy quark mass m^ s is then reconstructed as the 
invariant mass of these jets plus the opposite-sign lepton, of the three ones present. 
However, in the decay X — > Wt — > WWb, only half of the time the two jets and b 
quark will correspond to the top decay. We then set a cut 

140 GeV < m jlj2b < 210 GeV (32) 

to ensure that the event topology is consistent with the decay chain assumed. The 
number of events after the additional reconstruction conditions, including the cut in 
Eq. (1321) , can be found in Table [121 The m Y x distribution is shown in Fig. [TU For XX 
production we observe an endpoint around 500 GeV, which is not present for the other 
signals nor the SM background. Hence, if a signal is observed, template methods may 
be used to measure m x in this decay. Notice that a similar procedure to reconstruct the 
mass in BB decays is more difficult due to combinatorics, because the two W bosons 
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Figure 11: Visible reconstructed mass m^ s distribution of one of the heavy quarks (see 
the text). The luminosity is 30 fb -1 . 



from the heavy quark decay have opposite sign. Finally, it is worth remarking that 
the size of the signal itself would already give a strong hint that XX or BB pairs are 
produced, as it is apparent from the comparison of the numbers of events in Table [T2l 



5.3 Summary 

In this section it has been shown that the trilepton final state has very good sensitivity 
to TT, BB and XX production. Pair production of T, B and X quarks gives final 
states with three leptons with branching ratios not too small, and trilepton backgrounds 
can be significantly reduced with not very strong selection criteria, which keep most of 
the signals. The discovery potential found is similar to the dilepton final states, which 
have larger branching ratios but also larger backgrounds, but worse than in the single 
lepton one. However, the main interest of the trilepton final state is not heavy quark 
discovery but model discrimination, via the observation or not of several quark decays 
in this unique channel. Signals in the trilepton final state involve either the leptonic 
decay of a Z and a W boson, or of three W bosons. Hence, we have split the sample 
into two subsamples, one in which a Z candidate can be found (labelled as S Z') and 
the other one in which no such candidates can be identified (labelled as 'no Z') 

The £ ± £ ± £ =F (Z) final state is very interesting because it is very sensitive to TT and 
BB production, in the decay channels of Eqs. ([26]) . T and B singlets with a mass of 500 
GeV can both be discovered with a luminosity of 3.4 fb -1 . Discovery of a (T B) doublet 
requires 1.1 fb -1 (0.73 fb -1 ) in scenario 1 (2), and a (XT) doublet 0.72 fb -1 being the 

4 Notice that the transverse mass is not useful in this case because there are two more neutrinos 
from the other heavy quark decay. 
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main contribution to the signal from the T quark in the latter two cases. But this broad 
sensitivity brings an additional difficulty for the discrimination of signals: the three 
signal channels in Eqs. ( l26l) are interesting and it is necessary to identify event by event 
to which one it corresponds. In this sense, this final state is also very adequate because 
the kinematics of the decay chain can be fully reconstructed as there is only one light 
neutrino. With this purpose we have developed a likelihood analysis to discriminate 
among the three signal channels and then reconstruct the events accordingly. We 
have shown that the contributions of the three cascade decays in Eqs. ( 1261 ) can be 
disentangled, invariant mass peaks can be reconstructed, and the decays T — > Zt, 
T — > W + b, B — > Zb, B — > W~t can be established if sufficient luminosity is collected. 
This clean channel, in which the combinatorics is moderate, is also a good candidate to 
determine the chiralities of T, B quarks with the analysis of angular distributions in the 
semileptonic decay of the top quarks produced in T — ► Zt, B — > W~t. For example, a 
simple analysis presented shows that for 30 fb -1 the differences in a forward-backward 
asymmetry between a T quark singlet and a (T B) doublet in scenario 1 would amount 
to 2.4cx. In T — > Ht decays (seen for instance in the single lepton channel) the statistics 
is larger but the top quark polarisation is smaller and the reconstruction less clean than 
here. 

The final state without Z candidates is also very interesting because of its excellent 
sensitivity to XX and BB production, very similar to the one in the like-sign dilepton 
channel: 3.5 fb -1 for B singlets and 0.25 fb -1 for both (XT) doublets and the TB d2 
model, with the main contribution resulting from XX and BB production, respectively. 
Although the masses cannot be fully reconstructed in this final state, in XX production 
the X mass can be still determined from the endpoint of an invariant mass distribution. 
The presence of a signal for B quarks gives indirect evidence for the B — > W~t decay, 
which is absent for the B quark in a (BY) doublet. 

To conclude this section, it is worth mentioning that the trilepton final state is 
also very sensitive to heavy Dirac or Majorana neutrinos in singlet, doublet or triplet 
representations [35]. Those signals can be distinguished from heavy quark production 
because in that case the heavy neutrino can be observed as a peak in the invariant mass 
distribution of two opposite-charge leptons plus missing energy, with an additional 
peak in the distribution of the remaining lepton plus two jets. Scalar triplets also give 
trilepton signals (see for example Ref. [33]) but the like-sign dilepton invariant mass 
distribution displays a very sharp peak in scalar triplet production, which is of course 
absent in the case of heavy quarks. 
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6 Final state « 



This conspicuous signal can be produced in decays of B and X quark pairs when two 
same-sign W bosons decay leptonically, 

BB -> W~t W + i -> W~W + b W + W~b W ± -> ^i/, VT T -> gg' , 

xx -> iy + tiy-t -> -> ^i/, -> , (33) 

and also in decays of T quark pairs involving Z bosons 

Tf -> Zt W - & -> ZW^+fe^-fe Z -> £+r, WW -> £z/gg (34) 

when the opposite-charge lepton from the Z decay is missed by the detector (the charge 
conjugate channel is also included). Like-sign dilepton signals are relatively clean, their 
largest SM background being tinj in the semileptonic channel, where one of the two 
like-sign leptons results from a b quark decay. A very large source of £ ± £ ± events but 
for low lepton transverse momenta is bbnj, with a cross section of 1.4 ub (for a detailed 
discussion of like-sign dilepton backgrounds see Ref. [54]). For example, requiring only 
Pt > 15 GeV for the charged leptons the number of like-sign dilepton events from tinj, 
bbnj is around 25000 and 150000, respectively, for a luminosity of 30 fb _1 [54]. In 
order to reduce such backgrounds, we demand for event pre-selection (i) the presence 
of two like-sign leptons with transverse momentum pt > 30 GeV; (ii) the absence of 
non-isolated muons. The first condition practically eliminates bbnj while the latter 
reduces WZnj, which gives this final state when the opposite-charge lepton from Z 
decay is missed by the detector. The number of signal and background events at pre- 
selection can be read in Table H4l We note that it is sometimes claimed in the literature 
(without actually providing a proof) that SM backgrounds with charged leptons from 
b decays, namely tinj and bbnj, can be removed or suppressed to negligible levels by 
isolation criteria. However, recent analyses for supersymmetry searches performed with 
a full detector simulation [52] arrive at the opposite conclusion. This can already be 
seen at the level of a fast detector simulation. We show in Fig. [12] the minimum AR 
distance between the two charged leptons and the closest jet, for various signals and 
the tinj background. We observe that these variables, which in general do not bring 
a tremendous improvement in the signal to background ratio, are even inadequate in 
this case when the signals have many hard jets. 

For the like-sign dilepton final state we first perform a "discovery" analysis with 
conditions aiming only to improve the signal significance by reducing the background. 
Then, we impose additional requirements (which reduce the signal statistical signifi- 
cance) to try to reconstruct the event kinematics and detect heavy quark mass peaks. 
These two analyses are presented in turn. 
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79.3 


65.9 


BB (B s ) 


291.5 


170.1 


137.9 


Tf (TB dl ) 


156.3 


90.4 


74.8 


BB (TB dl ) 


368.6 


223.2 


172.8 


Tf (TB d jXT d ) 


263.8 


177.8 


149.1 


BB (TB d2 ) 


1737.4 


1122.8 


890.0 


XX (XT d ) 


1684.7 


1138.6 


900.4 


BB (BY d ) 


15.8 


5.6 


45.0 










YY (BY d ) 


2.0 


0.7 


0.2 


tinj 


1413 


43 


23 


WWnj 


245 


7 





Wtinj 


184 


47 


34 


WZnj 


1056 


9 


1 


Ztinj 


28 


9 


5 


WWWnj 


110 


11 


3 



Table 14: Number of events at the pre-selection, selection and reconstruction levels in 
the £ ± £ ± final state for the signals and main backgrounds with a luminosity of 30 fb -1 . 





Figure 12: Lego-plot separation between the leptons and the closest jet for several 
signals and the tinj background, at the pre-selection level (these variables are not used 
for event selection). The luminosity is 30 fb -1 . 



6.1 Discovery potential 

To evaluate the discovery potential for heavy quark signals we require for event selection 
(i) the presence of at least six jets, b tagged or not, with pt > 20 GeV; (ii) transverse 
momentum px > 50 GeV for the leading charged lepton l\\ (iii) missing energy /t> 50 
GeV; (iv) transverse energy larger than 500 GeV. The kinematical distributions of these 
variables at pre-selection are presented in Fig. [131 We also show for completeness the 
separate multiplicity distributions of light and 6-tagged jets. Notice that the maximum 
in the transverse energy distribution for the signals indicates that one or more heavy 
particles with a total mass around 1 TeV is produced. The number of signal and 
background events after our selection cuts is listed in Table [EH We observe that tinj 
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Figure 13: Light and 6-tagged jet multiplicity, total jet multiplicity, total transverse 
energy, transverse momentum of the leading lepton and missing energy. The luminosity 
is 30 ftr 1 . 



still amounts to one third of the total like-sign dilepton background after being reduced 
by the cuts. The corresponding discovery luminosities for each model can be found in 
Table [TH1 also indicating whether a mass peak can be reconstructed (see the next 
subsection). It is noticeable that not only XX and BB decays give like-sign dileptons 
but also TT decays, although the discovery luminosity in the T singlet model is much 
larger than for the rest. Finally, we note that for the B singlets and , (T B), (XT) 
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doublets the signals are much larger than the background and thus the uncertainty in 
the latter is not crucial for the evaluation of the discovery potential. 





L 


Rec. 




L 


Rec. 


T s 


17 fb" 1 


no 


TB d2 


0.23 fb" 1 


no 


B s 


4.1 fb" 1 


no 


XT d 


0.23 fb" 1 


m x 


TB dl 


1.5 ftr 1 


no 


BY d 




no 



Table 15: Luminosity L required to have a 5a discovery in the £ ± £ ± final state. A dash 
indicates no signal or a luminosity larger than 100 fb -1 . We also indicate whether a 
mass peak can be reconstructed in this final state. 



6.2 Heavy quark reconstruction 

In XX production the invariant mass of the quark decaying hadronically can be re- 
constructed from its decay products: a b quark and four jets from W decays. In order 
to do so, we restrict ourselves to events with at least one 6-tagged jet and four light 
(non-tagged) jets. The number of signal and background events after these additional 
reconstruction criteria is given in Table [HI The reconstruction is performed as follows: 

1. A 6-tagged jet is selected among the ones present. 

2. The four highest px light jets are grouped in two pairs J1J2, J3J4 trying to recon- 
struct two W bosons, the first one from the top quark decay and the second one 
from X -> Wt. 

3. The b jet is associated to the first light jet pair jij 2 to reconstruct a top quark. 

4. Among all the possible choices for the b jet and light jet combinations, the one 
minimising the quantity 

{m hh - M w f (m j3j4 - M w ) 2 (m jlj2b - m t f 
°w a w t 

is chosen, taking aw = 10 GeV, a t = 14 GeV. 

The reconstructed heavy quark mass m x is then defined as the invariant mass of the 
6-tagged and four light jets. These distributions are presented in Fig.[TH The long tails 
in the distributions of W and top reconstructed masses are mainly caused by wrong 
assignments. In particular, when only one 6-tagged jet is present in the event, half of 
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the times it corresponds to the b quark from the other heavy quark X in which the W 
bosons decay leptonically. Still, the heavy quark mass peak is clearly observed without 
the need of quality cuts on W and top reconstructed masses (which of course sharpen 
the mx peak). The rest of signals and the SM background do not exhibit any resonant 
structure, which shows that the above procedure does not introduce any bias. Note 
that a similar mass reconstruction cannot be achieved for BB or TT decays, due to 
the missing neutrino from each of the heavy quark decays. 

Although the cross section and mass reconstruction are consistent with the produc- 
tion of a XX pair we still can ask ourselves to which extent we can conclude that XX 
pairs are produced and not other possibility consistent with charge conservation and 
giving the same final state. We can establish XX production in two steps. First, we 
ensure that the extra boson from the heavy quark decaying hadronically, reconstructed 
from two jets j 3 and j 4 , is a W boson and not a Z. It is not surprising that the J3J4 
invariant mass distribution in Fig. HH has a peak around M w , since this is imposed in 
the reconstruction procedure. The question is then what would happen if, instead of 
choosing the pair of jets which best reconstruct a W boson, we chose the pair which 
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best reconstruct a Z. The comparison between both situations is shown in Fig. [THl 
including the two signals in the (XT) doublet (XX and TT). We observe that if 
we select the pair of jets with rrij 3 j 4 closest to M z the distribution is slightly shifted 
but the peak is maintained at Mw, and the heavy quark reconstruction is unaffected. 
Then, up to a more detailed study with a full detector simulation to confirm these 
results, it seems that the identity of the gauge boson from the heavy quark decay can 
be established. This leaves us with two options for the heavy quark: a B (charge —1/3) 
or X (charge 5/3) quark. (The opposite charges for antiquarks are understood.) 

Then, we examine the transverse mass distribution of the rest of particles in the 
event, restricting ourselves to events with two b tags for simplicity and imposing the 
quality cuts 

40 GeV < m jlj2 < 120 GeV , 

125 GeV < m jlj2b < 225 GeV (36) 
on the reconstruction (see Fig. [H]). We define the transverse mass as in Ref. [27], 

< = (4 ft + ^) 2 -(pD 2 , (37) 

with (E^ b ) 2 = (Px b ) 2 + m abi Pt = $t an d the transverse momenta of different particles 
summed vectorially. This distribution, shown in Fig. [16] for the relevant signals, has 
an edge around mx for the (XT) doublet signal, showing that the like-sign charged 
leptons and the b quark result from the decay of a 500 GeV resonance. Then, charge 
conservation and the absence of significant additional jet activity (which could be 
identified as additional partons produced in the hard process) implies the possible 
charge assignments (Qh, Qi) = ±(5/3, —5/3), ±(5/3, —7/3), ±(1/3, —5/3) for the heavy 
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quark decaying hadronically and leptonically, respectively. Of these three possibilities, 
the only one consistent with a small mixing of the third SM generation and the new 
quarks with the first two SM generations (so that b quarks are produced in bb pairs) is 
the first one, corresponding to XX production. 




Figure 16: Transverse mass distribution of the two charged leptons, a b jet and the 
missing energy. The luminosity is 30 fb^ 1 . 



6.3 Summary 

We have shown in this section that the like-sign dilepton final state has an excellent 
discovery potential for XX production: a (XT) doublet of 500 GeV could be discov- 
ered with only 0.23 fb -1 . This discovery potential is only matched by the trilepton 
and single lepton final states (but in the latter the quark observed is the T partner). 
Moreover, a heavy quark mass peak can be found in the invariant mass distribution of 
a reconstructed top quark and two extra jets, resulting from the hadronic decay of a 
W boson. Despite the fact that heavy quark charges cannot be directly measured (un- 
less the b jet charge is measured, which is very difficult at LHC), the detailed analysis 
of the event kinematics can eventually establish that the signal corresponds to XX 
production if this indeed is the case. 

For BB production the signals are also interesting and the discovery potential is 
also very good: 4.1 fb -1 for B singlets and 1.1 fb -1 , 0.23 fb -1 for (T B) doublets in 
scenarios 1 and 2, respectively. For BB production the heavy quark mass peaks cannot 
be reconstructed because each heavy quark has among its decay products an invisible 
neutrino. But the presence of a signal distinguishes a B singlet or [T B) doublet, which 
have decays B — > W~t, from the B quark in a (B Y) doublet which does not. This final 
state has some sensitivity to TT production, with one heavy quark decaying T — > Zt, 



44 



Ref. [27] 



Our analysis 



XX (XT d ) 
BB (TB d2 ) 
Background 



440 



424 



23 



470.6 



470.7 



116 



Table 16: Number of signal and background events for 10 fb 1 using the selection 
criteria in Ref. [27]. 

Z — ► £ + £~ and one of these leptons missed by the detector. The discovery potential 
is worse than in other final states, however, and the heavy quark masses cannot be 
reconstructed. 

We note that our results exhibit some differences with respect to previous work [27], 
due to two different sources: 

• We include a fast detector simulation and pile-up, and in order to reduce SM 
backgrounds we must apply tighter event selection criteria. For example, the 
event selection in Ref. [27] demands five jets with pseudo-rapidity \r]\ < 5. This 
is not sensible in the presence of pile-up, so we restrict our analysis to the central 
region of the calorimeter, \rj\ < 2.5. In fact, we have applied the cuts in Ref. [27] 
to our simulation, obtaining a similar signal efficiency but a background 5 times 
larger (see Table [TBI) . In particular ttnj, not considered there, amounts to 52 
events, about one half of the total SM contribution. Therefore, reducing the 
background requires stronger selection criteria which obviously reduce the signal 
as well. 

• As already indicated in section [3J our 5a discovery criterion is that (i) the signal 
statistical significance (possibly evaluated with Poisson statistics) is larger than 
5cr; (ii) the number of events is larger than 10. This second condition, which 
determines the limits on the XX and BB signals for this final state, is not 
included in Ref. [27]. Then, even for equal number of signal and background 
events, their limits are better. 

Our mass reconstruction method is also more involved and adapted to the more realistic 
conditions and the higher jet multiplicities found in our analysis. 

We finally comment on some other models giving the same final state. Like-sign 
dilepton signals without significant missing energy are characteristic of the presence of 
a heavy Majorana neutrinos [55], either in singlet or triplet SU(2)l representations (for 
a detailed comparison see Ref. [35]). Those models can be easily distinguished from 
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heavy quark production because in the heavy neutrino case (i) the missing energy is 
very small; (ii) two heavy resonances can be reconstructed, each one consisting of a 
charged lepton and two jets. On the other hand, like-sign dileptons with large missing 
energy are characteristic of heavy Dirac neutrinos in triplet SU(2)l representations [35], 
but in this case a resonance can be reconstructed with one charged lepton and two 
jets. Scalar triplet production also gives this final state but with the like-sign dilepton 
invariant mass displaying a sharp peak at the doubly charged scalar mass (see for 
example Ref. [33]). 



7 Final state £ + £~ 

This final state has large SM backgrounds which make it more difficult to observe 
positive signals by simply counting events with few selection criteria, as it is possible in 
the cleaner final states, and demand either a signal reconstruction to observe invariant 
mass peaks or an efficient background reduction. We ask for pre-selection the presence 
of (i) two opposite-charged leptons with px > 30 GeV; (ii) two 6-tagged jets with 
Pt > 20 GeV. Dilepton signals result from many signal decay channels, for example 

TT ^ ZtW~b ^ ZW + bW~b Z -> £ + r,W -> qq , 

Tf ^ZtVi^ ZW + b VW-b Z -> £ + T ,W -»• qq\V -»• qq/uu , 

BB -> ZbW + i -> ZbW + W-b Z -> £ + T,W -> qq' , (38) 

involving a Z — * £ + £~ decay (here V = Z, H), or 

Tf^W + bW~b W^lv, 

TT — > W + bVi — > W + bVW~b W -> Iv, V -> gg/i/P , 
fifi -> -> W~ W^+fe VT+^^fe -> £z/, 21^ -> gg , 

XX -> W^+t -> W - 6 2W -> £z/, 2iy -> gg' , 

y?^iy-&iyi w^^, (39) 

with two leptonic W decays. As it has been done in other final states, we separate the 
sample into two ones, one for events with a Z candidate (when the two charged leptons 
have the same flavour and invariant mass \rri£+£- — Mz\ < 15 GeV), and the other 
one with the rest of events, which do not fulfill one of these conditions. Backgrounds 
are also separated by this division: the ones involving Z production like Znj and 
Zbbnj mainly contribute to the former while tinj contributes to the latter. The total 
number of signal and background events at pre-selection level in both samples is given 
in Table [TFJ and the dilepton mass distribution for the signals in Fig. \17l 
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Total 


Z 


no Z 




Total 


Z 


no Z 


Tf (T s ) 


715.9 


179.6 


536.3 


BB (B s ) 


819.8 


393.4 


426.4 


Tf (TB dl ) 


799.4 


174.1 


625.3 


BB (TB dl ) 


907.5 


388.1 


519.4 


Tf (TB d jXT d ) 


1007.7 


341.6 


666.1 


BB (TB d2 ) 


1105.4 


55.4 


1050.0 


XX (XT d ) 


1147.4 


60.7 


1086.7 


BB (BY d ) 


902.5 


780.3 


122.2 










YY (BY d ) 


570.4 


25.8 


544.6 


tinj 


68493 


7464 


61029 


Z*/l*nj 


5245 


4875 


370 


tw 


2135 


212 


1923 


Zbbnj 


10132 


9807 


325 


tibb 


347 


38 


309 


Zccnj 


931 


883 


48 


Wtinj 


63 


4 


59 


Ztinj 


106 


88 


18 



Table 17: Number of events in the £ + £ final state for the signals and main backgrounds 
with a luminosity of 30 fb -1 , at pre-selection level. 




Figure 17: £ + £ invariant mass distributions for the six models in the £ + £ final state. 
The luminosity is 30 fb -1 . 



7.1 Final state £+£~ (Z) 

In the sample with \rri£+£- — M z \ < 15 GeV we first perform a generic analysis sensitive 
to T and B quarks, to obtain the discovery potential in this sample. Then, we perform 
a specific one aiming to detect the decay B — ► Hb (and thus the Higgs boson) in 
BB — ► ZbHb decays. This final state is interesting for Higgs boson discovery in models 
with doublets (B Y) where the only decays of the B quark are B — > Zb, B — ► Hb and 
the decay B — ► Hb cannot be reconstructed in the single lepton channel. 



47 



7.1.1 Discovery potential 

Here we demand for event selection (i) at least four jets with px > 20 GeV; (ii) trans- 
verse momentum pt > 50 GeV for the leading charged lepton l\\ (iii) transverse energy 
H T > 500 GeV. The kinematical distributions of the three variables are presented in 
Fig. EES Most signal channels, in particular those in Eqs. (l38j ). have at least four jets at 




4(10 



2500 

Figure 18: Kinematical distributions of variables used in selection criteria for the 

(Z) final state: light jet multiplicity, transverse momentum of the leading lepton and 

total transverse energy. The luminosity is 30 fb _1 . 

the partonic level. One exception is BB — ► ZbVb with V — » qqjvv^ in which additional 
jets are only produced by radiation or fragmentation. Still, this signal is sizeable after 
the multiplicity cut. The reader may also notice that the background could be further 
reduced by requiring for example Ht > 1 TeV. However, it is not clear whether this 
would indeed improve the signal observability. If the signal cannot be seen as a clear 
peak (or bump) over the background lineshape, its observation requires a simple event 
counting in which, if the background is large as it is our case, the background normali- 
sation uncertainty plays an important role. On the other hand, if the signal displays a 
peak the background can be in principle normalised from off-peak measurements and 
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its uncertainty will be smaller. The selection cuts made here represent a (conservative) 
compromise between having a manageable background and still observe the signal peak 
structure. 

As it was done for the trilepton channel, we build here a likelihood function to dis- 
criminate among the three signal channels in Eqs. (l38l) . building probability functions 
for three signal classes: (a) TT -> ZtWb; (b) TT -> ZtVt; (c) BB -> ZbWt. We 
generate high-statistics samples different from the ones used for the final analysis. We 
choose not to include a separate class for the background, because that would strongly 
bias it towards signal-like distributions and jeopardise the observation of reconstructed 
peaks. At any rate, the discriminant analysis implemented here rejects a large fraction 
of the background (which is classified as TT-like) when we concentrate ourselves on the 
BB signal. To build the discriminant variables we use an approximate reconstruction 
of the two W bosons decaying hadronically, choosing among the light jets (up to a 
maximum of 6) the four ones which best reconstruct two W bosons. Then, we use the 
same variables as in the trilepton channel but this time with two hadronic W bosons 
Wi, W2, ordered by transverse momenta as well as the 6-tagged jets b\ , 62- The result- 
ing distributions are presented in Fig. [T9l It is seen that the discriminating power is 
practically the same as in the trilepton channel, as it can be better observed in Fig. [20] 
which represents the likelihood function evaluated on the three class samples, giving 
the probabilities P a , P&, P c that the events correspond to each class. Table [TBI shows 
the performance of the likelihood function on the reference samples. Events in a class 
x are correctly classified if P x > P y , P z , where y, z are the other classes. 

Class P a >P b ,P c P b >P a ,P c P c >P a ,P b 



(a) 


0.59 


0.25 


0.16 


(b) 


0.23 


0.62 


0.15 


(c) 


0.17 


0.18 


0.65 



Table 18: Performance of the likelihood function on the event reference samples: 
fractions of events in each sample and their classification. Events in a class x are 
correctly classified if P x > P y , P z , where y, z are the other classes. 

The event reconstruction proceeds in the same way as in the trilepton channel (see 
section l5.1.2[ ) but replacing the leptonic W boson by a second W decaying hadronically. 
We use all jets pairings with a maximum of 6 light jets to construct two W bosons, and 
for events in class (b) we require at least 8 jets (6-tagged or not), otherwise the events 
are rejected. The number of signal and background events after reconstruction cuts 
and their distribution in the three classes is given in Table [lJB Notice that the total 
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Figure 19: Kinematical variables used to classify the three heavy quark signals in the 
£ + £~ (Z) final state. 



number of events includes those which are later rejected in reconstruction. We also 
remark that this signal discrimination based on topology brings an important "cleaning" 
of the background for events classified as BB, as we already have anticipated. 

We show in Fig. [21] the two most interesting signal peaks, those involving the 
decays T x — > Zt, common to classes (a, b), and B x — > Zb in class (c). These peaks 
are less biased by the reconstruction process although we can notice that, for example, 
the B — >• Zb singlet signals misidentified and reconstructed as T — » Zt have a small 
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Figure 20: Probability distribution functions for events in the reference samples. 




1250 



Figure 21: Reconstructed heavy quark masses in the £ + £ (Z) final state. The lumi- 
nosity is 30 fb _1 . 



bump around 500 GeV. Notice that the Bi peaks in the Zb invariant mass distribution 
are very sharp and significative, resulting in an excellent discovery potential for B 
quarks. For the other heavy quarks T 2 , B 2 with full hadronic decay the distributions 
are more biased, and the heavy quark peaks are almost equally well reconstructed for 
the signals as for the SM background: among the many jet combinations it is always 
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Total 


(a) 


(b) 


(c) 


Tf (T s ) 


137.3 


50.4 


56.5 


21.7 


Tf (TB dl ) 


135.6 


48.6 


56.6 


21.3 


Tf (TB d2 /XT d ) 


297.9 


67.4 


165.0 


40.5 


XX (XT d ) 


45.0 


14.2 


18.6 


2.8 


BB (B s ) 


220.4 


65.4 


31.1 


113.5 


BB (TB dl ) 


218.3 


61.1 


28.9 


121.3 


BB (TB d2 ) 


38.8 


11.3 


15.3 


3.3 


BB (BY d ) 


372.2 


134.8 


45.9 


180.9 


YY (BY d ) 


5.3 


1.9 


0.3 


2.7 


ttnj 


450 


113 


129 


16 


tw 


9 


4 


1 





Z*/l*nj 


181 


32 


86 


12 


Zbbnj 


335 


109 


74 


51 


Zccnj 


61 


24 


14 


8 


Ztinj 


65 


15 


28 


8 



Table 19: Number of signal events in the £ + £ (Z) final state at the selection level 
assigned to each event class. The luminosity is 30 fb _1 . 

possible, especially when the background involves top quarks and W bosons, to find 
one which is kinematically similar to the signal. These distributions are uninteresting 
and are not presented for brevity. Analogously, we do not present nor perform quality 
cuts on reconstructed W and top masses, which are very similar for the signals and 
backgrounds. We estimate the signal significance by performing a cut around the mass 
peaks, 

400 GeV < m Tl , m Bl < 600 GeV , (40) 

giving the numbers of signal and background events in Table l20l for completeness. The 
discovery luminosities obtained summing all signal contributions within a given model 
and combining the significances for the T\ and B\ peaks are presented in Table ED We 
find that the amount of work necessary to build the likelihood function and discriminate 
the different signals pays off, and the discovery luminosities achieved are quite small 
in some cases. The excellent result obtained for BB production in a (BY) doublet, 
where the B quark only decays in B — > Zb, B — > Hb, deserves a special mention: in 
most final states examined up to now the discovery potential for this model was rather 
limited but the opposite-charge dilepton one constitutes a remarkable exception. The 
reconstructed peak in the Zb invariant mass distribution shows the presence of a heavy 
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Ti (a, b) 


Bi (c) 




Ti (a, b) 


Bi (c) 


Tf (T s ) 


77.5 


15.1 


BB (B s ) 


55.7 


93.6 


Tf (TB dl ) 


74.2 


14.1 


BB (TB dl ) 


52.9 


103.7 


Tf (TB d2 /XT d ) 


159.6 


27.6 


BB (TB d2 ) 


13.2 


2.2 


XX (XT d ) 


18.7 


1.3 


BB (BY d ) 


101.3 


148.2 








YY (BY d ) 


1.4 


1.7 


ttnj 


69 





Zbbnj 


78 


28 


tw 


1 





Zccnj 


15 


4 


Z*/l*nj 


22 


8 


Ztinj 


22 


5 



Table 20: Number of signal events in the £ + £ (Z) final state at the Ti, Bi heavy quark 
peaks. The luminosity is 30 fb -1 . 





L 


Rec. 




L 


Rec 


T s 


22 fb" 1 


tut 


TB d2 


4.4 fb" 1 


ttit 


B s 


4.5 fb" 1 


m B 


XT d 


4.4 fb" 1 




TB dl 


2.4 fb" 1 


m T , m B 


BY d 


1.8 fb" 1 


m B 



Table 21: Luminosity L required to have a 5cr discovery in the £ + £ (Z) final state. 
We also indicate whether a mass peak can be reconstructed in this final state. 

B quark with charge —1/3. The same can be said regarding the Zt distribution and 
T quarks, although in this case the background is much larger and the observation in 
the trilepton final state is easier and cleaner. 



7.1.2 Discovery of B — > Hb 

We now concentrate ourselves on the process BB — > ZbHb. As selection criteria we 
only ask (i) the presence of four 6-tagged jets with px > 20 GeV, which is sufficient 
to practically eliminate all backgrounds, and (ii) less than four light jets, to remove 
the overlap between this final state and the previous one@ We give in Table [22] the 
numbers of events at pre-selection and selection for all signals and backgrounds. This 
final state has an excellent discovery potential for the Higgs boson with a (B Y) doublet 

5 This is not strictly necessary as long as we do not intend to combine the statistical sensitivities 
of both samples, but we include it for simplicity and in order to be conservative. Dropping the 
requirement on light jets the signals in this sample are larger and the discovery luminosity for the 
(B Y) doublet is reduced about a factor of two. 
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Pre. 


Sel. 




Pre. 


Sel. 


Tf (T s ) 


179.6 


3.1 


BB (B s ) 


393.4 


12.9 


Tf (TB dl ) 


174.1 


4.1 


BB (TB dl ) 


388.1 


14.4 


Tf (TB d2 /XT d ) 


341.6 


7.7 


BB (TB d2 ) 


55.4 


0.4 


XX (XT d ) 


60.7 


0.2 


BB (BY d ) 


780.3 


38.8 








YY (BY d ) 


25.8 


0.0 


tinj 


7464 





Z*/l*nj 


4875 





tw 


212 





Zbbnj 


9807 


6 


tibb 


38 


4 


Zccnj 


883 





Wtinj 


4 





Ztinj 


88 






Table 22: Number of events in the £ + £ (Z) final state at pre-selection and with the 
selection requirement of four 6-tagged jets. The luminosity is 30 fb -1 . 



L 


Rec. 




L 


Rec. 


T s 


no 


TB d2 




no 




no 


XT d 




no 


TB dl 30 fb" 1 


no 


BY d 


9.2 fb" 1 


m B , M H 



Table 23: Luminosity L required to have a 5a discovery in the £ + £~ (Z) final state 
with four b tags. We also indicate whether a mass peak can be reconstructed in this 
final state. 

(see Table l23l. and moderate for the TB dl model. In our estimations for the sensitivity 
we include a 20% systematic uncertainty in the background when necessary. 

The Higgs boson mass can also be reconstructed when it results from a B decay, 
doing as follows. 

1. We select a b jet to be paired with the Z boson candidate and reconstruct the 
heavy quark B\\ the other heavy quark B 2 is reconstructed from the three re- 
maining b jets. 

2. The combination minimising the mass difference (m^ 2 — is chosen. 

3. Among the three b jets from B 2 , we choose the two ones with minimum invariant 
mass to be the ones corresponding to the Higgs decay. 

The resulting reconstructed Higgs mass is shown in Fig. [22] (left). As we already 
have mentioned, this final state is most interesting for the (B Y) doublet which has 
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Figure 22: Reconstructed Higgs boson mass in the 
tags. The luminosity is 30 fb _1 . 



(Z) final state with four b 



a large signal and in which the Higgs mas peak can be clearly reconstructed with 
sufficient luminosity. For T — > Ht decays a small signal could be seen with different 
selection criteria but we do not address this here, since T — > Ht signals are far more 
interesting in the single lepton channel. Finally, one may wonder whether the "Higgs" 
peak results from the presence of a resonance or if it is merely a kinematical effect. To 
investigate this, we can use a different reconstruction by selecting the two b jets with 
smallest p?. The resulting distribution, shown in Fig. [22] (right), also displays a peak 
at the same place although the combinatorial background is larger in this case. 

7.2 Final state £ + £~ (no Z) 

In this final state the signals involve two W boson decays from different heavy quarks 
in general, and hence the heavy mass peaks are difficult to reconstruct except for BB 
production. The detection of a signal must then rely on event counting, which requires 
an efficient background suppression. We perform here two analyses: first a generic one 
aiming to discover the new quark signals, and then a specific one to reconstruct the 
heavy B quark mass. This mass reconstruction is useful for the TB^ 2 model, where 
the B — > Zb decay does not take place. 

7.2.1 Discovery potential 

For event selection we demand (i) transverse momentum px > 100 GeV for the sub- 
leading jet (6-tagged or not); (ii) transverse energy Ht > 750 GeV; (iii) the invariant 
mass of the highest-pr b jet b\ and each of the two leptons must be larger than the 
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top mass, taken here as 175 GeV. The first two conditions reduce backgrounds in 
general, while the third one strongly suppresses tinj production, where the b quarks 
and charged leptons result from top decays. The kinematical distributions of these 
variables are presented in Fig. [23J The number of events after these cuts can be read 
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Figure 23: Kinematical distributions of variables used in selection criteria for the £ + £~ 
(no Z) final state: transverse momentum of the second highest-pr jet, total transverse 
energy and invariant masses of the two leptons and the leading b jet. The luminosity 
is 30 fb -1 . 



in Table EH where we also include for better comparison the numbers of events at 
pre-selection. We point out that among the opposite-sign BB events produced in the 
decay 

BB — > W~tW + i — > W~W + bW + W~b 2W -> tv, 2W -> qq' , (41) 

those surviving the and cuts mostly correspond to leptonic decay of the 
two opposite-sign W bosons produced in B — > W~t, B — > W + t decays. Then, a 
mass reconstruction is difficult with these event selection criteria, which are anyway 
very efficient to reduce the tt background and observe a heavy quark signal. The 
luminosities required for 5a discovery are given in Table [251 Since in this final state 
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Pre. 


Sel. 




Pre. 


Sel. 


Tf (T s ) 


536.3 


236.6 


BB (B s ) 


426.4 


170.1 


Tf (TB dl ) 


625.3 


274.4 


BB (TB dl ) 


519.4 


202.7 


Tf (TB d JXT d ) 


666.1 


175.6 


BB (TB d2 ) 


1050.0 


233.3 


XX (XT d ) 


1086.7 


240.5 


BB (BY d ) 


122.2 


89.1 








YY (BY d ) 


544.6 


359.1 


tinj 


61029 


80 


Z*/l*nj 


370 


1 


tw 


1923 


14 


Zbbnj 


325 


6 


tibb 


309 


22 


Zccnj 


48 


1 


Wtinj 


59 


3 


Ztinj 


18 


6 



Table 24: Number of events in the £ + £~ (no Z) final state (discovery analysis) for 
the signals and main backgrounds with a luminosity of 30 fb -1 , at pre-selection and 
selection level. 





L 


Rec. 




L 


Rec. 


T s 


2.7 fb" 1 


no 


TB d2 


1.1 fb" 1 


m B 


B s 


9.3 ft" 1 


m B 


XT d 


1.1 fb" 1 


no 


TB dl 


0.83 fb" 1 


m B 


BY d 


0.87 fb^ 1 


no 



Table 25: Luminosity L required to have a 5a discovery in the £ + £ (no Z) final state. 
We also indicate whether a mass peak can be reconstructed in this final state. 

the background is still relatively important, the uncertainty in its overall normalisation 
affects the significance of the signals. We then include a 20% systematic uncertainty 
in our estimations in order to be more realistic. We point out the excellent sensitivity 
of this final state to TT, XX and YY production. The latter is specially important, 
because Y pair production only gives signals in the opposite-sign dilepton and single 
lepton channels, and the discovery of charge —4/3 quarks must be done in one of them. 

7.2.2 Heavy quark reconstruction 

In order to reconstruct the heavy B masses we drop from the selection criteria the £b 
invariant mass requirements, to allow for top quark semileptonic decays in Eq. ( T4"ll . 
The selection criteria in this case are: (i) the presence of four jets with px > 20 GeV; 
(ii) transverse momentum > 100 GeV for the sub-leading jet (6-tagged or not); (iii) 
transverse energy Ht > 750 GeV. The number of events are given in Table [26l We 
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Pre. 


Sel. 


Rec. 




Pre. 


Sel. 


Rec. 


Tf (T s ) 


■536.3 


209.3 


21.1 


BB (B s ) 


426.4 


206.4 


24.6 


Tf (TB dl ) 


625.3 


248.7 


27.0 


BB (TB dl ) 


519.4 


249.4 


36.2 


Tf (TB d2 /XT d ) 


666.1 


395.7 


42.4 


BB (TB d2 ) 


1050.0 


623.6 


122.4 


XX (XT d ) 


1086.7 


661.1 


127.2 


BB (BY d ) 


122.2 


53.6 


4.6 










YY (BY d ) 


544.6 


118.7 


11.0 


tinj 


61029 


1419 


139 


Z*/l*nj 


370 


2 





tw 


1923 


18 





Zbbnj 


325 


4 





tibb 


309 


28 


2 


Zccnj 


48 








Wtinj 


59 


6 





Ztinj 


18 


8 


2 



Table 26: Number of events in the £ + £~ (no Z) final state (reconstruction analysis) for 
the signals and main backgrounds with a luminosity of 30 fb _1 , at pre-selection and 
selection level, and including reconstructed mass cuts. 

notice that the tinj background is much larger here than in the previous discovery 
analysis (Table l24l. 

For events in which the W + W~ pair corresponds to the same heavy quark, the 
invariant mass of the W bosons decaying hadronically plus one of the b quarks will 
peak at rriB- The heavy quark mass reconstruction is done analogously as for the X 
quark in the like-sign dilepton channel in section 16.21 and the reconstructed heavy 
quark mass tub is then defined as the invariant mass of the b quark and four light jets 
selected. These distributions are presented in Fig. EH We point out that, in contrast 
with the X quark reconstruction in the like-sign dilepton channel, half of the events 
in BB decays have opposite-sign leptons resulting from different heavy quark decays. 
The heavy B quark mass peak is reasonably well reconstructed as it is shown in the 
last plot of Fig. EH but the distribution is quite similar for a X quark. Then, although 
the presence of a signal would be apparent, the observation of a clear peak and the 
discrimination among these two possiblities is rather difficult, even more in the presence 
of a large ti background^ We then apply quality cuts to improve the reconstruction 
and reduce the background, 

60 GeV < m jlj2 < 100 GeV , 
60 GeV < m j:ij4 < 100 GeV , 
125 GeV < m jlhb < 225 GeV . (42) 

6 For higher masses the background suppression is more efficient via transverse energy requirments, 
and the B quark peak may be easier to reconstruct [56]. 
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Figure 24: Reconstructed masses of the two W bosons, the top and the heavy quark. 
The luminosity is 30 fb _1 . 



The number of events after these cuts is given in Table [26l With these cuts, the 
reconstructed mass for the TB d2 model (displaying a peak at m B ) and the (XT) 
doublet (without peak) are quite different, as it can be seen in Fig. [25] (left). Both 
possibilities could also be distinguished in the presence of background, as shown in 
the right panel. The reconstruction of a peak in the Wt invariant mass distribution 
shows that the heavy quark has charge —1/3 or 5/3. Unfortunately the transverse 
mass of the B quark with leptonic W decays does not display a clear endpoint, and the 
direct identification of the quark charge is not possible as for X quarks in the like-sign 
dilepton channel. 



7.3 Summary 

Despite the a priori large opposite-sign dilepton backgrounds, this final state turns 
out to have an excellent sensitivity to heavy quark signals. These signals and the SM 
backgrounds result from either the leptonic decay of a Z or of two W bosons. Therefore, 
as it was done in other final states, it is advantageous to divide this final state into two 
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Figure 25: Left: reconstructed mass of the heavy quark for the two largest signals. 
Right: the same, including the SM background. The luminosity is 30 fb -1 . 



subsamples, with or without a Z candidate. 

The dilepton sample with a Z candidate has a huge background from Znj produc- 
tion which can be practically removed by asking for the presence of two 6-tagged jets 
and four light jets. Once that SM backgrounds are manageable, we have implemented a 
likelihood function which discriminates among the three processes in Eqs. ([3*51) . classify- 
ing events with a good efficiency. The event reconstruction is then performed according 
to this classification. For events identified as BB — > ZbWt, the Zb invariant mass dis- 
tribution displays a sharp peak which signals the presence of a charge —1/3 quark. For 
TT production a peak can be observed in the Zt invariant mass distribution but the 
background is larger in this case. The discovery luminosities are small, as it can be 
seen in Table EU The case of B quarks in the (BY) doublet (1.8 fb -1 for 500 GeV 
quarks) deserves a special mention, since these quarks are much harder to see in other 
final states. For this model, the observation of the decay B — > Hb, with H — > bb is 
also possible if we concentrate on final states with four b tags, and 5cr discovery could 
be possible with a luminosity below 10 fb -1 for mg = 500 GeV and Mh = 115 GeV. 
This decay is interesting not only to establish the identity of the B quark but because 
it is a possible discovery channel for the Higgs boson if such doublets exist. 

In the subsample without Z candidates we have performed two different analyses, 
first a generic one which achieves the best signal significance and then a specific one 
to reconstruct the heavy B quark mass in B — > W~t decays. The background is again 
important but the largest one, tinj, can be practically removed by requiring invariant 
masses m$, > m t , so that the charged leptons and b quarks cannot result from a top 
quark decay (this requirement must be dropped in the reconstruction analysis). After 
background suppression, this final state also offers an excellent discovery potential 
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for the pair production of T, X and Y quarks. For example, 5<r significance can be 
achieved with a luminosity around 1 fb -1 for the four models with quark doublets. 
The sensitivity to Y quarks is especially important because they only have decays 
Y — > W~b, and their detection can only be performed in the opposite-sign dilepton 
and single lepton final states. 

Finally, we have performed the reconstruction of the B quark mass in B — > W~t — ► 
W~~W + b decays (or the charge conjugate), with both W bosons decaying hadronically. 
This is specially interesting for the TB d , 2 model where the B — ► Zb decay does not 
take place. In this analysis the > m t requirement must be dropped in order to 
keep the events which actually display a peak in the invariant mass of four light jets 
plus a 6-tagged jet. With adequate reconstruction quality cuts a clear peak could 
be observed distinguishing BB and XX production, but the quark charge cannot be 
directly measured. 

8 Final state £± 

Single lepton signals result from heavy quark pair decays when one of the W bosons 
(up to four are present, depending on the channel) decays leptonically and the rest of 
W, Z and Higgs bosons decay hadronically. Hence, single lepton signals benefit from a 
large branching ratio. This final state is fundamental to establish whether the T — > Wb 
decay takes place (this can also be seen in the trilepton final state but needs about ten 
times more luminosity, see section [5. 1.21) . Due to the large size of the signals, one can 
also look for subsamples with high b jet multiplicities to establish the decays T — > Ht, 
B — ► Hb. Without being completely exhaustive, we will perform three analyses in this 
section. The first one, in the single lepton channel with exactly two b tags, is devoted 
to the search of the T — > Wb decay. The second one, in a sample with four b tags, 
allows to search for T — > Ht and B — > Hb. The third one, requiring six b tags, is very 
useful to look for T — > Ht in the models where this decay has an enhanced branching 
ratio. For event pre-selection we require (i) one charged lepton with pt > 30 GeV; (ii) 
at least two b jets with pr > 20 GeV; (iii) at least two light jets also with px > 20 GeV. 
The total number of events for the signals and main backgrounds at pre-selection are 
given in Table EZJ as well as the numbers in the three subsamples. 
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Total 


2b 


46 


66 


Tf (T s ) 


9415.3 


5797.8 


874.7 


26.5 


Tf (TB dl ) 


10064.4 


6172.9 


931.8 


29.0 


Tf (TB d jXT d ) 


11782.9 


5294.3 


1873.4 


112.1 


XX (XT d ) 


9213.8 


7506.3 


172.9 


0.4 


BB (B s ) 


6535.6 


4460.1 


412.9 


2.5 


BB (TB dl ) 


7021.4 


4802.3 


434.2 


2.8 


BB (TB d2 ) 


9193.4 


7484.4 


164.6 


0.0 


BB (BY d ) 


2146.1 


1399.6 


150.7 


0.7 


YY (BY d ) 


7444.1 


6588.5 


58.7 


0.3 


tinj 


965051 


902205 


1629 





tw 


31920 


30280 


38 





ttbb 


4355 


2287 


423 


2 


titi 


27 


12 


1 





Wnj 


38185 


37236 


14 





Wbbnj 


20634 


19920 


16 





Wtinj 


654 


592 








Z/ynj 


3397 


3314 








Zbbnj 


4874 


4715 


5 






Table 27: Number of signal and background events in the £ ± final state at the pre- 
selection level, and in the four subsamples studied. The luminosity is 30 fb -1 . 

8.1 Final state £ ± (26) 

In this channel, the processes we are mainly interested in are 

Tf -> W + b W~b WW -> £uqq' , 

YY -> W~b W + b WW -> iuqq' . (43) 

There are other decay channels for T quarks, involving T — > Zt and T — > Ht, which 
give the same final state, but they are suppressed by the selection criteria. Some 
production and decay channels of B and X quarks also give single lepton signals, for 
example 

BB -> W't W + i -> W~W + b W + W-b 3W -> qq\ 1 W -> tv , 

XX -> W + tW~t -> W + W + bW~W-b 3W -> gg', 1W -> £z/ , (44) 

but they are rather difficult to see due to their greater similarity with ti production. 
Then, we concentrate our analysis on the channels in Eqs. f T43l) . The selection criteria 
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applied for this are: (i) transverse momentum pr > 150 GeV for both b jets; (ii) 
transverse energy H T > 750 GeV; (iii) the invariant mass of both b jets and the charged 
lepton must be larger than the top mass, taken here as 175 GeV. The distributions of 
the relevant variables are shown in Fig. [261 The first condition is inspired by the specific 
decays in Eqs. ( 1431 ). (The dependence on the pt cut is not very strong, and we have 
choosen 150 GeV for simplicity.) The transverse energy requirement is a general one 
to look for high mass states, not very optimised for these heavy quark masses. The 
invariant mass requirements are extremely useful to reduce the tinj background (still 
some events remain due to mistag of the charm or light jets), and allows to improve our 
results over previous analyses for TT production in this channel [24,25]. However, it 
reduces the BB and XX signals for which some of the decay channels have the charged 
lepton and b quarks both resulting from a top quark. We give in Table [28] the number 
of events at selection, also including the ones at pre-selection for better comparison. 
We observe the excellent background reduction achieved with these simple selection 
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Figure 26: Kinematical distributions of variables used in selection criteria for the £ ± 
(26) final state: transverse momentum of the subleading b jet, total transverse energy 
and invariant masses of the charged lepton and the two b jets. The luminosity is 30 
fb- 1 . 
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Pre. 


Sel. 


Peak 




Pre. 


Sel. 


Peak 


TT (T s ) 


5797.8 


796.6 


538.8 


BB (B s ) 


4460.1 


391.7 


217.7 


TT (TB dl ) 


6172.9 


795.0 


551.7 


BB (TB dl ) 


4802.3 


352.8 


193.4 


TT (TB d2 /XT d ) 


5294.3 


123.8 


55.2 


BB (TB d2 ) 


7484.4 


186.5 


84.7 


XX (XT d ) 


7506.3 


165.5 


74.8 


BB (BY d ) 


1399.6 


307.0 


198.8 










YY (-Bid) 


6588.5 


1974.3 


1508.8 


tinj 


902205 


299 


117 


Wbbnj 


19920 


125 


52 


tw 


30280 


68 


34 


Wtinj 


592 


12 


1 


ttbb 


2287 


15 


6 


Z/'ynj 


3314 


2 


1 


tut 


12 


1 





Zbbnj 


4715 


24 


9 


Wnj 


37236 


49 


22 











Table 28: Number of signal and background events in the (26) final state at the 
pre-selection and selection level, and at the reconstructed mass peak. The luminosity 
is 30 fb^ 1 . 

criteria, especially with the £b invariant mass cuts: the tinj background is reduced by 
a factor of 3000, while the YY signal is kept at one third. The TT signal is reduced to 
one seventh because there are contributing channels other than TT — > W + bW~b, and 
those are quite suppressed by the event selection necessary to reduce tinj. 

The TT and YY signals are reconstructed by choosing the best pairing between b 
jets and reconstructed W bosons: 

1. The hadronic W is obtained with the two jets (among the three ones with largest 
Pt) having an invariant mass closest to M w . 

2. The leptonic W is obtained from the charged lepton and the missing energy with 
the usual method, keeping both solutions for the neutrino momentum. 

3. The two heavy quarks Q = T,Y are reconstructed with one of the W bosons 
and one of the b jets. We label them as Qi^, corresponding to the hadronic and 
leptonic W, respectively. 

4. The combination minimising 

(mj% - M w f (mj% - M w f (mjg - m^) 2 
°W °W °Q 

is selected, with aw = 10 GeV, <jq = 20 GeV. 
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Figure 27: Reconstructed heavy quark masses at the selection level in the £ ± (26) final 
state. The luminosity is 30 fb _1 . 



We present the reconstructed mass distributions at the selection level in Fig. [271 With 
the criteria applied the peaks are very good even for the TT signals which involve 
several competing decay chains giving up to six b quarks. In particular, the pt cut on 
the subleading b jet suppresses the T — > Ht decays selecting only T — >■ Wb, in which 
we are interested. The signal significance can be estimated by performing the invariant 
mass cuts 

350 GeV < m 1 ^ 2 < 650 GeV . (46) 

The numbers of signal and background events at the peak, defined by the above mass 
windows, can be found in Table EH The high signal significance achieved for the case of 
T and Y quarks implies an excellent discovery potential, summarised in Table [291 We 
include a 20% systematic uncertainty in the estimations in all cases. For T quarks the 
discovery luminosities are rather small except for the TB d2 and (X T) doublet model, 
where the T — > Wb decay does not take place. For YY production the discovery po- 
tential is even better, because the Y — > Wb channel is the only one present. For B and 
X quarks the signals are smaller, of the size of the background itself, implying a signal 
significance below 5cr even for large luminosities due to the background normalisation 
uncertainty assumed. However, in case that a signal is detected in other final states it 
should be possible to detect also BB and XX signals in the single lepton channel with 
two b tags using a dedicated and optimised analysis. The heavy quark reconstruction 
as a peak in Wb invariant mass distributions implies that it has either charge 2/3 or 
—4/3. Both possibilities cannot be distinguished unless the b jet charge is measured, 
which is very difficult. However, for the models considered in this paper a strong hint 
is offered by the signal size itself, which is much larger for YY production than for XT, 
and the observation of T — > Zt in the opposite-sign dilepton and trilepton final states 
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L 


Rec. 




L 


Rec. 


T s 


1.1 fb" 1 


rriT 


TB d2 


- 


no 


B s 




no 


XT d 




no 


TB dl 


0.60 fb" 1 


rriT 


BY A 


0.18 fb" 1 


m Y 



Table 29: Luminosity L required to have a 5a discovery in the £ ± (26) final state. A 
dash indicates no signal or a luminosity larger than 100 fb -1 . We also indicate whether 
a mass peak can be reconstructed in this final state. 

establishes the quark charge. 

Finally, it is worth mentioning that T quarks singlets and those in the TB Al model 
(which have the same decay channels) could in principle be distinguished by the W 
helicity fractions [57], but for large tut the W bosons in T — > W + b are mainly longitu- 
dinal and the difference between a left- and right-handed WTb coupling is washed out. 
For a 500 GeV T singlet we have F L ~ 0.05, F ~ 0.95, F R ~ 0, while for the T quark 
in a (TB) doublet F L ~ 0, F ~ 0.95, F R ~ 0.05. With 500 events for a 30 fb" 1 , the 
statistical error ~ 1/y/N in angular asymmetries, etc. is expected to be around 5%, 
of the order of the difference between the two models. Systematic uncertainties in the 
measurement of helicity fractions can also be important [53]. 

8.2 Final state £ ± (46) 

The main interest of this final state, in addition to the discovery of new quarks, lies 
in the observation of the decays T — * Ht, B — > Hb, which would also allow an early 
discovery of a light Higgs boson if new vector-like quarks exist [26,30,32, 58]. The 
relevant decay channels are 

iuqq' , 

iuqq', V — * qqj vv , 
ivqq . (47) 

The first and last channels in the above equation give exactly four b quarks in the 
final state, while the second gives up to six b quarks. For TT production alone, it 
has been shown [26] that the discrimination between the first two decay chains is very 
involved in final states with only four 6-tagged jets, because the signals are actually 
not very different. The situation is worsened if additional B quarks exist, for example 
in models introducing a (T B) doublet. Here we implement a discriminating method 



Tf -> Ht W~b -> HW + bW-b H -> 66, WW 

TT -> HtVi^ HW + b VW~b H -> 66, WW 

BB -> HbW + i -> HbW + W'b H->bb,WW 
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Figure 28: Total transverse energy distribution for the signals and backgrounds in the 
for the tr (46) final state. The luminosity is 30 fb _1 . 

based on a likelihood analysis similar to the ones used in previous sections. In this 
final state, however, discrimination is less efficient than for multi-lepton signals due to 
the combinatorics resulting from presence of four 6-tagged jets and their association 
to the other particles present in the event. We use high statistics reference samples 
for three event classes (a,b,c) corresponding to the three decay channels in Eqs. fT4"7"l) . 
As selection criteria for this analysis we demand (i) four 6-tagged jets with px > 20 
GeV; (ii) transverse energy Ht > 750 GeV. The distribution of this variable for the 
signals and the SM background can be found in Fig. EH The variables used in the 
discrimination are obtained using a preliminary reconstruction of a top quark and the 
Higgs boson: 

1. The hadronic W is obtained with the two jets (among the three ones with largest 
Pt) having an invariant mass closest to M w . 

2. The leptonic W is obtained from the charged lepton and the missing energy with 
the usual method, keeping both solutions for the neutrino momentum. 

3. The top quark is reconstructed with one of the W bosons and one of the four b 
jets, selecting the ones which give an invariant mass closest to the nominal top 
mass. The W boson and b quark selected are labelled as W 2 and b 2 . 

4. The Higgs boson "candidate" is obtained from the two b jets, among the three 
remaining ones, which have the minimum invariant mass. 

5. The remaining W boson and b jet are labelled as W±, b\. 
The interesting variables for signal discrimination are: 
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• The light jet multiplicity. 

• The Wibi invariant mass , which peaks around my in class (a). 

• The Hbi invariant mass, which peaks at mg in class (c). 

• The invariant mass. This corresponds to m T in class (a), but the distri- 
bution is very similar for the other decay channels. 

• The WiW 2 b2 invariant mass, which is m# in class (c) but does not differ much 
for events in the other decay channels. 

The normalised distributions of these variables for the reference event samples are 
presented in Fig. EHJ together with the resulting probability distributions P a ,b,c that 
the events belong to a given class. Comparing with the trilepton and dilepton final 
states, we find that the separation among channels is much less clean. This is reflected 
in Table [3(2 which collects the fractions of events correctly and incorrectly classified. 

Class P a >P b ,P c P b >P a ,P c P c >P a ,Pb 



(a) 


0.53 


0.17 


0.30 


(b) 


0.28 


0.47 


0.25 


(c) 


0.27 


0.15 


0.58 



Table 30: Performance of the likelihood function on the £ (46) event reference samples: 
fractions of events in each sample and their classification. Events in a class x are 
correctly classified if P x > P y , P Zl where y, z are the other classes. 

We give in Table [31] the numbers of signal and background events at the selection 
level, and their classification (as in the previous sections, we select the class which has 
the highest probability). An important remark here is the presence of XX and YY 
signals, as well as BB ones in the TB d2 model, which involve only two b quarks at 
the partonic level. Events with four 6-tagged jets result from the mistag of the charm 
quarks from W decays, and light quarks to a lesser extent. The presence and size 
of these signals illustrates the relative importance of mistags in the processes we are 
interested in: for TT and BB events with four b jets it may well happen that only 
three of them correspond to true b quarks and one is a charm quark from a W decay. 
This, added to the kinematical similarity of the signals and the several possibilities in 
b jet assignments, makes the separation among the channels difficult. The luminosities 
required for 5a discovery are collected in Table [32], summing all contributions in a 
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Figure 29: Kinematical variables used to classify the three heavy quark signals in the 
£ ± (46) final state, and the resulting probability distributions for events in the reference 
samples. 
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Total 


(a) 


(b) 


(c) 


Tf (T s ) 


836.8 


342.5 


260.4 


233.9 


Tf (TB dl ) 


886.5 


363.9 


286.4 


236.2 


Tf (TB d2 /XT d ) 


1780.7 


509.9 


841.9 


428.9 


XX (XT d ) 


167.3 


44.9 


86.5 


35.9 


BB (B s ) 


396.8 


119.6 


64.3 


212.9 


BB (TB dl ) 


416.4 


119.7 


67.2 


229.5 


BB (TB d2 ) 


160.0 


43.0 


83.1 


33.9 


BB (BY d ) 


146.1 


62.0 


10.5 


73.6 


YY (-Bid) 


57.9 


28.1 


4.9 


24.9 


tinj 


404 


122 


228 


54 


tw 


5 


3 





2 


ttbb 


158 


47 


66 


45 


titi 


1 








1 


Wnj 


1 





1 





Wbbnj 


3 





1 


2 


Zbbnj 


1 








1 



Table 31: Number of signal and background events in the £ ± (46) final state at the 
selection level assigned to each event class. The luminosity is 30 fb -1 . 





L 


Rec. 




L 


Rec 


T s 


0.70 fb" 1 


to t , M h 


TB d2 


0.16 fb" 1 


no 


B s 


1.9 fb" 1 


m B , M H 


XT d 


0.16 fb" 1 


no 


TB dl 


0.25 fb" 1 


to t , m B , M H 


BY d 


6.2 fb" 1 


no 



Table 32: Luminosity L required to have a ha discovery in the £ ± (4b) final state. We 
also indicate whether a mass peak can be reconstructed in this final state. 

given model and combining the significance of the three classes (a,b,c). A systematic 
uncertainty of 20% is included in the estimations. 

Comparing with Ref. [26], we observe that the discovery luminosity is significantly 
smaller than the values quoted there. The reasons for this difference are: (i) for 
consistency with the other channels we are giving here the statistical significance of 
the signal (including the Higgs boson) compared to the "only background" hypothesis, 
while in Ref. [26] we compared the "Higgs" and "no Higgs" hypotheses in the presence 
of new quarks, for which the significance is lower; (ii) multi-jet SM backgrounds in 
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Ref. [26] were quite pessimistically overestimated, and with our new evaluation with 
updated tools and improved matrix element-parton shower matching they turn out to 
be smaller; (iii) the new likelihood classification performed here, with the subsequent 
statistical combination of channels, also improves the significance. 

We finally address the heavy quark and Higgs boson reconstruction, which depends 
on the decay channel in which events are classified. We only perform the reconstruction 
of events classified in the first and third classes, because those in the second class have 
six b jets at the partonic level. 

Class (a): TT — > HtWb — > HWbWb. Events identified as resulting from this decay 
chain are reconstructed using this procedure: 

1. Two light jets are selected to form the hadronic W, labelled as Wh- If there are 
only two light jets these are automatically chosen; if there are more than two, 
only up to three (ordered by decreasing p?) are considered. 

2. The leptonic W (labelled as Wl) is obtained from the charged lepton t and the 
missing energy, in the way explained in previous sections. Both solutions for the 
neutrino momentum are kept, and the one giving best reconstructed masses is 
selected. 

3. Two b jets are selected among the ones present, to be paired with Wh and Wl, 
respectively. 

4. The top quark is reconstructed from one of the Wb pairs, and its parent heavy 
quark Ti from the top quark and the two remaining b jets. 

5. The other heavy quark T 2 is reconstructed from the remaining Wb pair. 

6. Among all choices for b and light jets and all possible pairings, the combination 
minimising the quantity 



is selected, with a t = 14 GeV, gt = 20 GeV. After the final choice, the Higgs is 
reconstructed from the two b jets not assigned to the W bosons. 

Class (c): BB — > HbWt — > HbWWb. The reconstruction of this channel proceeds 
through the same steps 1 — 2 as in the previous two channels, and then: 

3. One b jet is selected and paired with one of the two W bosons to form a top 
quark, and with the other W to form the heavy quark B 2 . 
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(48) 



71 




Figure 30: Reconstructed heavy quark and Higgs masses in the (46) final state. 



4. The three remaining b jets then reconstruct the heavy quark B\. 

5. The combination minimising 

(mj% - M w f (mj% - My/) 2 « - m t ) 2 (mg - mg 2 c ) 2 

^ + ~2 + ~2 + ~2 ' i 49 ) 

with a B = 20 GeV, is finally selected. Among the three b jets corresponding to 
Bi , the two with the minimum invariant mass are chosen to reconstruct the Higgs 
boson. 
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The results are presented in Fig. [301 For brevity we do not include the reconstructed 
W boson and top quark distributions, which have good peaks by construction. We ob- 
serve that for TT signals in class (a) the T\ >2 peaks are well reconstructed, as the Higgs 
boson peak. The same happens for BB signals in class (c): the B\ y2 peaks are clear and 
the Higgs boson peak is observable. However, the intriguing fact is that, for TT and 
BB signals included in the "wrong" class (respectively, (c) and (a)) the reconstruction 
procedure produces peaks which are as sharp as those for the signals "correctly" clas- 
sified, except for the Higgs peaks. This point deserves a detailed discussion. Clearly, 
if an event is assigned to a given decay channel in Eqs. ( jjTj) based on its likelihood, it 
is because its kinematics is quite compatible with that decay channel. Then, it is not 
so surprising that, for example, if a BB event is classified as TT based on its topol- 
ogy, when it is reconstructed as TT it rather looks as a TT event. For the reader's 
illustration we present in Fig. [3D several distributions for the production of TT and 
BB events in the case of singlets in all decay channels. In the upper part we show 
the normalised W±bi and W 2 b 2 H invariant mass distributions assigned to class (a) For 
BB events (incorrectly classified) the distributions display peaks very similar to the 
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Figure 31: Comparison between kinematical distributions for correctly and wrongly 
classified events (see the text). 
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ones for TT events correctly included in this class. In the lower part of this figure we 
plot the WiW 2 b 2 and Hbi distributions for TT and BB events in class (c). The peaks 
are quite similar for TT (wrong classification) and BB (correct). Therefore, we can 
conclude that distinguishing TT and BB signals in this channel is a more demanding 
task, and the multi-leptonic channels are much more appropriate for that. Fortunately, 
all these difficulties in signal discrimination do not affect the discovery potential, which 
is excellent for this final state. 

8.3 Final state t± (66) 

The single lepton final state with six 6 jets allows a clean reconstruction of the decay 

TT -> HtHi^ HW + b HW'b H -> 66, WW -> £uqq' , (50) 

with H — > 66, which seems impossible if only four jets are tagged. This final state 
is most interesting for the models in which the decay T — > Ht is enhanced and the 
66 signal is larger. We do not impose any further selection criteria apart from having 
six 6-tagged jets with pt > 20 GeV, which defines the sample studied. The number 
of background events is given in Table [331 This final state is extremely clean, and 





Sel. 




Sel. 


TT (T s ) 


26.5 


BB (B s ) 


2.5 


TT (TB dl ) 


29.0 


BB (TB dl ) 


2.8 


TT (TB d2 /XT d ) 


112.1 


BB (TB d2 ) 


0.0 


XX (XT d ) 


0.4 


BB (BY d ) 


0.7 






YY (BY d ) 


0.3 


tibb 


2 







Table 33: Number of events in the £ ± (66) final state at selection level. The luminosity 
is 30 fb" 1 . 

discovery could be made merely by an event counting. The 5cr discovery potential 
for the different models is given in Table [34J summing all signal contributions. The 
discovery potential (in models with T quarks, when a signal is produced) is determined 
by the requirement of having at least 10 signal events. The background normalisation 
in this case has little effect on the significance, because for the discovery luminosities 
it is rather small. 

The event reconstruction can be easily done despite the large combinatorics from 
the six 6 jets. The procedure is similar to the ones used in other final states: 
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L Rec. 

T s 11 fb" 1 m T , M H 

B s no 

TB dl 9.4 fb" 1 m T , M H 



L Rec. 

TB d2 2.7 fb" 1 m T) M H 

XT d 2.7 fb" 1 m T , M H 

BY A - no 



Table 34: Luminosity L required to have a 5a discovery in the £ ± (66) final state. A 
dash indicates no signal or a luminosity larger than 100 fb -1 . We also indicate whether 
a mass peak can be reconstructed in this final state. 



1. Two light jets (among the three ones with largest px) are selected to form the 
hadronic W, labelled as Wh- 

2. The leptonic W, labelled as Wl, is obtained from the charged lepton £ and the 
missing energy. 

3. Two b jets are selected among the ones present, to be paired with the two W 
bosons to reconstruct the top quarks decaying hadronically and semileptonically 
(t H and t L ). 

4. The four remaining b jets are grouped in pairs to reconstruct the two Higgs 
bosons, Hi and H 2 . 

5. The two heavy quarks T\ (corresponding to Wh) and T 2 (with Wl) are recon- 
structed from a top quark plus a Higgs boson. 

6. Among all choices for b and light jets and all possible pairings, the combination 
minimising the quantity 

{m^ H - M w f (mj% - M W ) 2 (mg - m t ) 2 {m\f - m t ) 2 

°W a W t t 

°\ o\ [ ' 

is selected, with on = 20 GeV. 



The results are presented in Fig. [32l We omit for brevity the W boson and top quark 
reconstructed masses, which good have peaks at My/ and m t by construction. It is 
seen that for the TB d2 and (X T) models the heavy quark and Higgs peaks are quite 
good and, with a moderate luminosity, they would give evidence for the T — > Ht decay 
(and, in particular, for the production of a Higgs boson). For T singlets and the TB dl 
model the signals are smaller and they would require more luminosity, not only to be 
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Figure 32: Reconstructed heavy quark and Higgs masses in the ^ (6b) final state. 



discovered but also to reconstruct the peaks. We point out that an important difference 
with the 46 final state is that six b jets can only be produced (up to mistags of charm 
quarks) from the decay in Eq. (1501 ). Hence, the model identification is cleaner here. 
We also note that for TT production within the TB d2 and (X T) models the signal in 
this final state is almost four times larger than in the £ + £~ (Z) one with four b tags, 
so this final state is best suited to detect T — > Ht. 



8.4 Summary 

The single lepton final state offers the best heavy quark discovery potential for all the 
models studied, due to the large signal branching ratios. To achieve this result, an 
efficient reduction of the large backgrounds from ttnj, Wnj and Wbbnj is necessary. 
We have concentrated on three different subsamples in this final state, with exactly 
two, four and six b jets. Clearly, the discovery potential will improve further including 
the samples with three and five b jets, which have not been considered here for brevity. 

The final state with two 6-tagged jets is the best suited for the discovery of YY 
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production, which only requires 0.18 fb for a 500 GeV quark. It is also very good for 
T singlets (1.1 fb -1 for the same mass) and (T B) doublets in scenario 1 (0.60 fb -1 ). In 
this final state the T and Y masses can be reconstructed as peaks in the Wb invariant 
mass distributions. The identity of the quarks cannot be established unless the b jet 
charge is measured: the decays T — > W + b and Y — > W + b both give a W + boson 
plus a 6-tagged jet, and the jet charge measurement is necessary to discriminate both 
possibilities. It is interesting to point out that all kinematical distributions are the 
same for T and Y quarks, including various angular asymmetries which can be built in 
the W and top quark rest frames (see for example Ref. [59]). The only discrimination 
between both possibilities comes either indirectly, from the cross section measurement 
(about three times larger for YY in this final state, after including efficiencies and cuts) 
or directly, via the observation of T — > Zt in the dilepton or trilepton final states. 

The final state with four b jets is the best one for the discovery of T quarks in 
either of the models considered. For heavy quark masses of 500 GeV, 5a discovery 
only requires 0.7 fb -1 for singlets, 0.25 and 0.16 fb -1 for the TB dl and TB d2 models, 
respectively, and 0.16 fb -1 if the T quark belongs to a (XT) doublet. For B singlets 
the discovery potential is also the best one, with 1.9 fb -1 . (For these results we have as- 
sumed a light Higgs, as suggested by precise electroweak data, taking M H = 115 GeV.) 
This process could also be a discovery channel for the Higgs boson in the presence of T 
or B quarks. We have gone beyond the signal observation and studied the discrimina- 
tion among TT and BB signals in this final state, which makes sense because they are 
both present in general for the case of the (T B) doublet. The separation is very diffi- 
cult, for several reasons: (i) the combinatorics from the presence of four b jets; (ii) two 
different decay chains contribute in the case of TT; (iii) the signals are kinematically 
not very different; (iv) the possibility of charm quark mistags. We have implemented 
a likelihood method to separate TT and BB signals, which has a reasonable efficiency 
if we bear in mind all these difficulties. After the TT and BB events are classified, the 
kinematics can be reconstructed according to the decay channel expected in each case, 
and sharp peaks are obtained in all cases, although the rate of "wrong" classifications 
is sizeable and a better discrimination between T and B quarks can be achieved in the 
trilepton final state. 

It is also interesting to remark the excellent discovery potential for (X T) doublets 
in this final state: only 0.16 fb -1 for heavy quark masses of 500 GeV. The discovery 
potential in this channel is similar but better than in the like-sign dilepton and trilepton 
channels (0.23 and 0.25 fb -1 , respectively) altough in those final states the main signal 
contribution comes from the charge 5/3 quark X and here it is the T quark which gives 
most of the four-6 signal. The same results are obtained for the TB d2 model which 
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includes a T quark with the same decay modes: 5a discovery of T is possible with 
0.16 fb -1 in the single lepton final state with four b tags, while the B quark can be 
discovered in the like-sign dilepton and trilepton final states with 0.23 and 0.25 fb -1 , 
respectively. 

Finally, the sample with six b jets has also been studied. In this final state the 
decay TT — > Ht Hi can be cleanly determined and peaks reconstructed without con- 
tamination of other decay modes. The signals are small, however, except for the TB d2 
and (XT) models, for which T — > Wb does not take place and thus T — > Ht has a 
larger branching ratio. The discovery potential is rather good for these models, 2.7 
iV 1 for m T = 500 GeV. 

9 The roadmap to top partner identification 

We summarise in Table [351 the discovery luminosities for the six models in the different 
final states examined. The comparison among them clearly shows that the single lepton 
channel (either with two or four b jets) offers the best discovery potential for new quarks. 
In the case of doublets, the signals may correspond to one or both members, as it is 
explained in detail in the summary at the end of each section, where it is also indicated 
whether heavy quark masses can be reconstructed. We now discuss case by case how 
the different models would be discovered and identified. 





T s 




B s 




TB dl 


TB d2 




XT A 




BY d 




£+£+£-£- (ZZ) 






24 ftr 
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23 fo" 
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10 fb" 
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£+£+£-£- (Z) 
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50 fb" 
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-i 
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Table 35: Luminosity required to have a 5a discovery in all final states studied. 
A T singlet or the T quark in the TB dl model: They would be discovered in the 
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single lepton final state with two or four 6 jets. In the 26 final state, the peaks in the 
Wb invariant mass distributions give evidence of the charged current decay but do not 
identify completely the new quark (it could be a charge —4/3 quark Y, although the 
cross section would not be consistent with that hypothesis). In the 46 final state, the 
peak in the Ht distribution, with H reconstructed from two b jets exhibiting a peak at 
M H , gives quite strong hints of the T — > Ht decay, although B quarks also give signals 
not very different. (The 6b final state does not have this ambiguity but the observation 
of T — > Ht requires much larger luminosity.) The best confirmation of its nature 
comes with a little more luminosity in the £ ± £ ± £ =F (Z) final state, with the observation 
of a peak in the Zt invariant mass distribution. This peak also establishes that the 
quark has charge 2/3. The analysis of the charged lepton distribution in the top quark 
rest frame for the subset of events in which the top quark decays semileptonically 
can discriminate between a T singlet with a left-handed WTb coupling and the T 
quark in a (T B) doublet for which the coupling is right-handed. With 30 fb _1 the 
differences found in the forward-backward asymmetry would amount to 2.4<r, and a 
better sensitivity is expected by using a more sophisticated analysis with a fit to the 
complete distribution. 

A T quark in the TB d2 model or in a (XT) doublet: in the single lepton final 
state with two b quarks it does not exhibit peaks in the Wb invariant mass distribution 
because the decay T — > W + b does not take place and probably the signal is very difficult 
to separate from the tinj background. In the 4b sample, however, the signal is very 
large and clean, and the quark is seen in the decay T — > Ht. With small luminosity, 
a signal should be also visible in the 66 sample. This quark also has enhanced decays 
T — > Zt, from which the quark charge is determined, and the signals in the trilepton 
and opposite-sign dilepton final states with a Z candidate are 2 — 3 times larger than 
expected for a T singlet. 

A B singlet or the B quark in the TB dl model: they would be discovered in the single 
lepton final state with four 6 jets. However, its discrimination from a T quark might 
not be very clear due to combinatorics and the signal similarities. With practically the 
same luminosity, the B quark would appear as a sharp peak in a Zb invariant mass 
distribution in the i^i^i* (Z) final state. This would determine the quark charge, 
and would be confirmed by an opposite-sign dilepton signal. The evidence for the 
B — > W~t decay comes from the same trilepton final state, and the charged lepton 
distribution in the top semileptonic decays would in principle probe the chirality of the 
WtB coupling, but the statistics is smaller than for T quarks. Indirectly, evidence for 
the B — > W~t decay results from the presence of £ ± £ ± .£ =F (no Z) and ^ ± £ ± signals, also 
observable with small luminosity. 
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A B quark in the TB d2 model: it would be discovered in the like-sign dilepton and 
£ ± £ ± £ =F (no Z) final states with similar luminosities. An indirect indication of the quark 
identity, in comparison with a charge 5/3 quark X, would be given by the absence of 
the reconstructed mass peaks and endpoints which are present for XX production. 
Signals in the £ + £~ (no Z) final state are also interesting, not only because of the 
good discovery potential but also because the mass reconstruction is possible in the 
hadronic decay B — > W~t — > W~W + b (or the charge conjugate) with a moderate 
luminosity. This mass reconstruction is important because in this model the B quark 
does not have decays B — > Zb and thus trilepton and opposite-sign dilepton signals 
with a Z candidate are absent (see the paragraph above). Single lepton signals with 
four 6-tagged jets, which are very significant for other models, are also absent for the 
B quark in this model but can be produced by its T partner and are kinematically not 
very different. 

A B quark in a (B Y) doublet: it does not give trilepton signals without a Z candi- 
date nor like-sign dilepton ones. On the other hand, it gives large opposite-sign dilepton 
signals with a Z candidate with a sharp peak in the Zb invariant mass distribution, 
from which the quark charge is determined. With five times more luminosity, this is 
also done in the four lepton final state with two Z candidates. The decay B — > Hb can 
be seen in the £ + £~ (Z, 46) final state, also with larger luminosity. 

A charge 5/3 quark X: it would be simultaneously discovered in the like-sign dilep- 
ton and £ ± £ ± ^ =F (no Z) final states. In the former, the invariant mass can be recon- 
structed and the quark identity (i.e. that it has charge 5/3) can be established under 
reasonable assumptions. The mass could also be determined from the trilepton final 
state with an endpoint analysis to confirm the quark identity. A signal due to this 
quark should also be visible in the four lepton final state without Z candidates. 

A charge —4/3 quark Y: it would be discovered in the single lepton final state with 
two b jets. The peaks in the Wb invariant mass distributions would give evidence of 
the charged current decay and indirect evidence of its nature: the signal is three times 
larger than for a T singlet, for example. A clean signal in the £ + £~ (no Z) final state 
would also be visible, larger than for a T quark. On the other hand, all the signals 
characteristic of a T — > Zt decay would be absent, in the trilepton and dilepton final 
states with a Z candidate, for example. 

In more complicated scenarios with several singlets or doublets the signals would 
add up but it still would be possible to identify the new quarks with a thorough 
examination of all final states. For example, with a T singlet plus a (XT) doublet in 
which both charge 2/3 quarks have similar masses (which cannot be experimentally 
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resolved), the T — > W + b decay would be seen in the single lepton final state, and the 
dilepton and trilepton signals involving T — > Zt would be much larger than the ones 
corresponding to just one T quark. Of course, the X member of the doublet would also 
be detected. The same argument applies to a combination of a B singlet and a (B Y) 
doublet. The simplest discrimination between T, B singlets and the T, B quarks in the 
TB dl model is by the presence of two partners almost degenerate in mass. Still, one 
may imagine a situation in which a T and a B singlet almost degenerate were present. 
This scenario would be distinguished from the TB dl model with the analysis of angular 
distributions, especially for T — > Zt decays. 

The discrimination of vector-like singlets and doublets from a fourth sequential 
generation with new quarks t', b' is also easy, because the latter would give different 
signatures (see Refs. [6, 60] for recent reviews). Oblique corrections prefer a mass 
difference m t > > my ~ 60 GeV [5], so the t' quark would decay either t' — > W + b (as 
a T or Y) or t' — > W + b' if the mass difference is larger. In the first case, the absence 
of t' — > Zt would distinguish it from a T singlet, and the absence of a degenerate B 
partner from a (B Y) doublet. In the second case, the decay between new heavy quarks 
would prove the non-singlet nature of both. Regarding the b' quark, for my > m t + M w 
the decay b' — > W~t would dominate, distinguishing this quark from a B singlet in 
either model. The discovery potential for fourth generation quarks [61,62] is similar to 
the models studied here. 

To summarise, the analyses carried out in this paper show that the single lepton 
final state offers the best discovery potential, and is the one in which new vector-like 
quark signals would be first seen. Searches in the dilepton and trilepton channels would 
soon confirm a possible discovery, and with a luminosity around five times larger all 
the decay modes of the new quarks would be observed in these channels, establishing 
the nature of the new quarks. In some models four lepton signals could be sizeable 
and detectable as well and, in any case, these should be investigated as a further test 
of the models. 

10 Conclusions 

In this work we have investigated in detail the LHC discovery potential for pair pro- 
duction of new vector-like quarks in five models: T or B singlets of charge 2/3, —1/3 
respectively, and (T B), (XT), (BY) doublets of hypercharge 1/6, 7/6, —5/6, restrict- 
ing ourselves to the case that new quarks mainly couple to the third generation, as it 
is expected from the SM quark mass hierarchy. In the case of (T B) doublets we have 
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distinguished two scenarios: that both heavy quarks have similar mixing with the top 
and bottom quark (model TB dl ) and that the mixing of the top with its heavy partner 
is much larger than for the bottom quark (model TB d2 ), as expected from the mass 
hierarchy m t m& and from indirect precision data. Using a dedicated Monte Carlo 
generator Protos [43] we have computed all signal contrubutions involving all heavy 
quark, gauge and Higgs boson decay channels. With a fast detector simulation of sig- 
nals and backgrounds we have examined twelve final states which would give evidence 
of the presence of new quarks, with one to four charged leptons in different kinematical 
regions and several b jet multiplicities. 

We have identified the final state with one charged lepton plus two or four b jets 
as the most sensitive one for new quark searches. Nevertheless, model discrimination 
requires the observation or exclusion of the different heavy quark decay channels. To 
achieve this goal, the dilepton and trilepton final states are essential. These final states 
have also good sensitivity to heavy quark signals, and with a luminosity at most five 
times larger than in the single lepton channel the 5a observation would be possible and 
the heavy quarks might be identified. The reconstruction of mass peaks would also be 
possible when a sufficient number of events is collected. In our simulations we have 
taken heavy quark masses of 500 GeV, focusing on early discoveries at LHC. We have 
obtained an excellent discovery potential for all models: 0.70 and 1.9 fb -1 for T and 
B singlets, respectively; 0.25 and 0.16 fb -1 for the TB dl and TB d2 models; and 0.16, 
0.18 fb -1 for the (XT) and (BY) doublets. It is also interesting to know the mass 
reach for higher integrated luminosities. With a simple rescaling it can be seen that in 
the single lepton channel alone and a luminosity of 100 fb -1 heavy T, B singlets with 
masses up to 800 and 720 GeV respectively can be discovered with 5a significance, 
while for the doublets the reach is higher: 850 GeV and 900 GeV for the (T B) doublet 
in the two scenarios considered, 900 GeV for (XT) and 820 GeV for (BY). For higher 
masses the experimental detection of heavy quarks can also be done using jet mass 
measurements [56] but model discrimination would follow similar strategies as outlined 
here. 

We have also obtained an excellent potential for the discovery of the new quarks in 
decay channels containing a Higgs boson, especially in the final state with one charged 
lepton and four 6-tagged jets. For heavy quark masses of 500 GeV, the discovery 
luminosities are 0.16 fb -1 for the TB d2 and (XT) models, 0.25 fb -1 for TB dl and 
0.70, 1.9 fb -1 for T and B singlets, respectively. These luminosities are much smaller 
than the ones required for a light Higgs discovery in the SM. Indeed, it is well known 
since some time [32, 58] that vector-like quark production can be a copious source of 
Higgs bosons and, if such quarks exist and the Higgs boson is light, its discovery would 
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possibly happen in one of these channels. For a heavier Higgs with different decay 
modes the analyses presented here (relying on the leading decay H — > bb) must be 
modified accordingly. Nevertheless, the determination of the other modes like T — > 
W + b, T — > Zt, etc. would still be done in the same way as presented here, with few 
modifications. 

In the summaries given at the end of each section we have compared the multi- 
lepton signals produced by new quarks with those arising from heavy leptons. Both 
possibilities for new fermions are easily distinguished by the different reconstructed 
mass peaks and the common presence of b jets for quarks, which in lepton pair produc- 
tion only result from H — > 56, Z — > bb decays. Interestingly, a more general difference 
among models introducing new quarks and leptons is that the latter give signals which 
are more "multi-leptonic": for heavy leptons the trilepton signatures are usually the 
ones with the highest significance, while for heavy quarks the single lepton one is the 
most sensitive. This is not unexpected, since heavy lepton decays give SM leptons plus 
a gauge or Higgs boson, while heavy quarks give SM quarks instead. In the minimal 
supersymmetric standard model where squark and gluino pair production is large, it 
is also found [52] that, although multi-lepton signatures are important, the final state 
with best discovery potential is the one with a charged lepton or large missing energy 
plus jets. 

Finally, it is worth pointing out that, although in this work we have restricted 
ourselves to heavy quark pair production, electroweak single production can also have 
a large cross section depending on the mass and couplings of the new quarks. These 
interesting processes are the only ones in which the heavy quark mixing with the SM 
sector can be measured because the pair production cross section is determined by the 
heavy quark mass alone and the heavy quark total width is likely to be very difficult 
to measure. Further model discrimination is also possible in single production, in 
particular from the study of cross sections and angular asymmetries, and it will be 
addressed elsewhere. 
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A Feynman rules 



We give in Tables [3614401 the Feynman rules used for our matrix element evaluations for 
heavy quark pair production TT, BB, XX, YY and single production Tbj, Bbj, Xij, 
Ttj, Ybj. The rest of vertices are the same as in the SM, including QCD interactions 
of the new quarks. 




T t t T 



Table 36: Feynman rules for T singlet electroweak and scalar interactions with the 
third generation. 



w + 




B b b B 



Table 37: Feynman rules for B singlet electroweak and scalar interactions with the 
third generation. 
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Table 38: Feynman rules for (T B) doublet electroweak and scalar interactions with 
the third generation. 




Table 39: Feynman rules for (XT) doublet electroweak and scalar interactions with 
the third generation. 
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w + n 



i ^ ir X&g (tr-bPl + m b P R ) 

ZM\y 



B b 
H 




-iTT-XhB-fPR 
2cw 



i 77^— X *bB( m bPL + m B P R ) 



B 



Table 40: Feynman rules for (B Y) doublet electroweak and scalar interactions with 
the third generation. 
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